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Abstract. Gyrokinctic theory is based on an asymptotic expansion in the small 
parameter e, defined as the ratio of the gyroradius and the characteristic length of 
variation of the magnetic field. In this article, this ordering is strictly implemented 
to compute the electrostatic gyrokinetic phase-space Lagrangian in general magnetic 
geometry to order e 2 . In particular, a new expression for the complete second- 
order gyrokinetic Hamiltonian is provided, showing that in a rigorous treatment 
of gyrokinetic theory magnetic geometry and turbulence cannot be dealt with 
independently. The new phase-space gyrokinetic Lagrangian gives a Vlasov equation 
accurate to order e 2 and a Poisson equation accurate to order e. The final expressions 
are explicit and can be implemented into any simulation without further computations. 
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1. Introduction 

Gyrokinetics [T] has proven a very useful tool to study turbulence in the core of fusion 
devices, making kinetic simulations of turbulent fluctuations possible in reasonable 
computational times [2j El HJ 13 El [7]. Its main advantage is averaging over the 
gyrofrequency time scale without losing the effect of the finite size of the gyroradius that 
is of the order of the typical wavelength of the turbulence. To perform this average, it is 
necessary to assume certain orderings that in the electrostatic limit can be summarized 
as 

B(r) with V ~ — 
L 

1 Id 
(plT.t) with V _l ~ — , b ■ V ~ — , — ~w 

TK ' ; p L dt 

~6<1, (1) 



Q L Mv 2 t 



where f(r,t) is the electrostatic potential, B(r) is the magnetic field, u is the 
characteristic frequency of the turbulent fluctuations, L is a characteristic macroscopic 
scale, v t , p = v t /Q and Q = ZeB/Mc are the thermal speed, the gyroradius and the 
gyrofrequency of the species of interest, Ze and M are the charge and the mass, and 
e and c are the magnitude of the electron charge and the speed of light. Since this 
article is about electrostatic gyrokinetics, we have assumed that the magnetic field is 
stationary and its characteristic length of variation is of the order of the macroscopic 
length L. The ordering in ([T]) implies that the electrostatic potential fluctuates with 
some characteristic frequency u and has a strong gradient perpendicular to the magnetic 
field, on the order of the inverse of the gyroradius, whereas its gradient parallel to the 
direction of the magnetic field, b = B/fi, is on the order of the inverse of the larger scale 
L. The frequency u of the turbulence is usually much smaller than the gyrofrequency, 
making the gyrokinetic average over the gyromotion valid. We have employed the small 
parameter e ~ u/H < 1 to make this explicit. In most fusion experiments, the ratio 
of the gyroradius and the macroscopic length is another small parameter that we also 
order as e. It is easy to see that in drift wave turbulence, for which the characteristic 
frequency is u ~ v t /L, the quantities u/fl and p/L are indeed of the same order. More 
importantly, to obtain the typical gyrokinetic formalism, it is necessary to order the 
electrostatic potential as small compared with the characteristic energy of the particles. 
This assumption is necessary to prove that the gyromotion of the particles is circular 
to lowest order. The most common gyrokinetic ordering assumes that the parameter 
Zeip/Mvt is comparable to e, as is done in (JT|). In this way, magnetic geometry effects 
such as the V-B and curvature drifts, of order p/L, are allowed to be comparable to 
the turbulent E x B drift, of order Ze^f/Mv^. Ordering these effects so that they are 
comparable is very important in, for example, the core of tokamaks, where the curvature 
of the magnetic field lines is believed to be the most important drive for the turbulence 
[8]. The ordering in ([[]) contains the simplest assumptions that are still interesting, 
but it can be extended to include components of the potential that have perpendicular 
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gradients of the order of the inverse of the macroscopic length L [91 [10]. The results 
that we present in this article can be easily extended to some of these more general 
orderings, but we leave this for future work. 

There are different techniques to obtain gyrokinetics (and for that matter, drift 
kinetics [HI [12j [13], of which gyrokinetics is a natural extension). On the one hand, 
it is possible to obtain the gyrokinetic equation by working iteratively on the Vlasov 
equation [HI U5J [TBI ELT]- We will call these iterative methods. On the other hand, 
it is possible to use phase-space Lagrangian/Hamiltonian methods that solve order by 
order for the motion of the particle in a given electromagnetic field, uncoupling the 
gyromotion from the slower time scales [181 H3 EH EU [22], [231 EH [25]. Once the motion 
of the particle is known, the Vlasov equation is simply obtained by its characteristics. 
We will call these Lagrangian methods. Both procedures are asymptotic expansions 
in the parameter e, and give equivalent equations order by order, but the Lagrangian 
methods have the advantage of giving the equations in a form that exactly conserves 
some energy-like quantities. This property may be very important for the global, full / 
simulations that are being developed [251 EZ1 ESI El]- To have an energy-like invariant 
and at the same time obtain equations of motion and equations for the electromagnetic 
fields that are the same to first order in e as those obtained with the iterative procedure, 
it is necessary to carry the expansion in e to higher order. For example, in a slab [23], it is 
necessary to obtain the Hamiltonian to second order in e. The second order piece of the 
Hamiltonian, quadratic in the electrostatic potential <p, gives second-order corrections 
to the equations of motion and hence it is in principle negligible to first order. However, 
if the lowest-order quasineutrality equation that contains a linear term in ip is employed, 
the second-order correction to the Hamiltonian must be kept to obtain an energy-like 
invariant. 

The complete calculation to order e 2 has not been done for a general static magnetic 
field in either formalism so fajj]. In the most common Lagrangian formulation [25], the 
calculation is done in two steps: first, the turbulent electromagnetic fields are ignored 
and only the background magnetic field is considered, giving the drift kinetic equation; 
in the second step, the turbulent electromagnetic fields are added and the corresponding 
corrections are calculated. Consider the case in which the magnetic field does not vary 
in time, i.e. electrostatic gyrokinetics. In the first step, the equations are expanded 
in the small parameter e ~ p/L, whereas in the second step, they are expanded in 
e v ~ Zeip/Mvf. The expansion in e is only performed to first order because the next 
order results are very tedious to calculate. The expansion in e v is continued to second 
order because the pieces quadratic in (p are needed to have an energy- like invariant. In 
the expansion in e v , the fact that there has been a previous expansion in e is ignored, 
and as a result the terms of order ee v are never calculated. The missing terms of order 
e 2 and ee^ are comparable to the terms of order e 2 according to the gyrokinetic ordering 
in (1), making this expansion consistent only when ^> e. In addition, since the cross- 

| In the particular case of a constant magnetic field the calculation to order e 2 was given in |23j . 
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terms that contain both the background magnetic field and the turbulent electrostatic 
potential, of order ee v , are always neglected when the two-step method is presented, 
it is not obvious how to calculate them following that procedure. In this article, we 
present the complete phase-space Lagrangian calculation with the standard gyrokinetic 
ordering (1), emphasizing the self-consistent calculation of the terms of order e 2 and ee^. 
In the gyrokinetic equations that result from the new Lagrangian, the magnetic geometry 
effects and the fluctuating potential appear together in the second-order terms, showing 
that geometry and turbulence cannot be separated and dealt with independently. Our 
main result is the explicit expression for the second-order gyrokinetic Hamiltonian given 
in equations (11321) . (11331) . (I134p and (11351) . It clearly exhibits the interplay between 
geometry and turbulence inherent to gyrokinetic theory, possessing terms of three types: 
terms quadratic in the electrostatic potential, terms that include both the electrostatic 
potential and the magnetic geometry, and terms that are purely geometrical. 

At this point, it is fair to wonder about the motivations beyond formal coherence 
to carry out the expansion consistently to second order in e. Keeping the second order 
piece of the Hamiltonian that is quadratic in the electrostatic potential is necessary for 
the conservation of an energy-like invariant, as already noted above. When the other 
second-order terms computed in this article are included, they have two effects: (i) the 
gyrokinetic Poisson's equation is modified by the effect of the non-uniform magnetic field 
on the gyro-orbits, and (ii) the equations of motion are modified to second order. Both 
of these effects are not conventionally kept in gyrokinetic formulations, but they may be 
crucial for conservation of momentum. Conservation of momentum in full / gyrokinetic 
formulations has been the center of a recent controversy [T7J [301 ED E2J E21 El]- By 
assuming a gyroBohm level of turbulent transport of momentum at long wavelengths, 
Catto and one of us, FIP, have argued that to recover with a full / model the correct 
transport of toroidal angular momentum in a tokamak, it is necessary to have gyrokinetic 
Fokker-Planck and Poisson's equations correct to third order in e in the high flow 
ordering, for which the average velocity of the ions Vi is of the order of the ion thermal 
speed Vti, and correct to fourth order in the low flow ordering, for which Vi ~ evu <C v t i- 
In the case of slab gyrokinetics, a consistent calculation of the transport of momentum 
in the low flow ordering requires the third order Hamiltonian [33, 34J. The requirements 
for a system with general geometry are still to be sorted out, and to do so it is necessary 
to study the new terms presented here and terms of even higher order. 

In addition to the issues raised for full / simulations, the formulation presented here 
will be very useful for 5f approaches to momentum transport in tokamaks in the low 
flow ordering. Reference [35] presents a formulation of this problem in the electrostatic 
limit that requires the minimum number of modifications to existing Sf simulations. 
The most important conclusion in [35] is that the turbulent pieces of the distribution 
function and the electrostatic potential have to be calculated to an order higher in 
e than usual because the contribution to momentum transport from the lowest order 
pieces vanishes due to symmetry arguments [361 EZ]- These symmetry arguments do 
not hold if the higher order terms of the gyrokinetic equation are considered. To avoid 
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calculating most of the next order corrections to the gyrokinetic Vlasov and Poisson's 
equations, reference has to resort to a subsidiary expansion based on the fact that in 
many tokamaks the poloidal component of the magnetic field is much smaller than the 
toroidal component. The new contribution to the Hamiltonian that we calculate here 
gives the self-consistent higher order contributions to the gyrokinetic Vlasov equation 
for the first time. Only a higher order gyrokinetic Poisson's equation is then lacking to 
obtain a complete 5f formulation in the low flow ordering that does not require a small 
poloidal magnetic field; this will be the subject of a future publication. 

The rest of this article is organized as follows. In Section [2] we write the non- 
dimensional phase-space Lagrangian of a particle in an electromagnetic field. The 
normalization shows explicitly the standard gyrokinetic ordering ([T]). In the first part 
of Section [3] we review the phase-space Lagrangian approach to gyrokinetics to help 
the understanding of the calculation that follows. In the second half of this section we 
proceed to obtain the gyrokinetic Lagrangian to second order in our expansion parameter 
e. As mentioned above, this is our main result. In Section H] we obtain the Vlasov 
equation from this Lagrangian, and in Section \5\ we discuss the consequences of this 
formulation for Poisson's equation. The new Vlasov and Poisson's equations presented 
here are correct to second and first order in e, respectively. We should remark that 
in the limit where the electrostatic potential has a scale of variation much larger than 
the gyroradius of the species of interest, our gyrokinetic equations provide the highest 
order guiding-center equations that we are aware of. In Section O we borrow tools from 
classical field theory to obtain Poisson's equation in a different way. We prove that there 
is an energy-like invariant and we discuss the stringent conditions on the equations to 
actually conserve it in a simulation. We finish with a discussion of our results and the 
future lines of research in Section [7J The Appendices contain the most cumbersome 
parts of the calculation as well as some material included for completeness. Finally, we 
would like to stress that in this article we have given all our results in an explicit form 
that can be directly implemented in a computer code. 

2. Normalized Lagrangian 

The phase-space Lagrangian for the motion of a particle of mass M and charge Ze in 
an electromagnetic field is given by 



Here A is the vector potential that is defined such that B = V x A. Notice that the 
phase-space Lagrangian depends on the position of the particle r, its velocity v, the time 
derivatives of both the position and the velocity, r = dr/dt and v = dv/dt, and the 
time t. For convenience, we will sometimes denote the phase-space coordinates {r,v} 





with the Hamiltonian 




(3) 
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as {X a } & a=1 = X = {r,v}. We use the superscript x in the Lagrangian ([2]) because it 
is a function of the phase-space coordinates X. 

The equations of motion are obtained by finding the stationary points of the 
action cr x [r(t), v(t)] = dt £ x (r(t), v(t), r(t), v(t), t) with respect to variations of the 
functions r(t) and v(t) subject to the constraints r(t = t ) = r , v(t = t ) = v , 
r(t = t x ) = rx and v(t = ix) = V!. From this procedure we obtain six equations of 



motion, namely 
d 
dt 

and 



^(V r £ x ) = V r £ x (4) 



|(Vv£ x ) = V v £ x . (5) 

This differs from the standard Lagrangian formalism where the Lagrangian function 
depends only on r, r and t. Actually, the phase-space Lagrangian formalism can be 
viewed as a variational formulation of Hamilton equations (see, for example, [2S])- In 
plasma physics, it was first applied by Littlejohn to guiding-center dynamics in [2T] . 

The Lagrangian ([2]) is non-dimensionalized using the characteristic thermal velocity 
of the species of interest Vt, the characteristic length L^ 1 ~ |V(ln|A|)| and the 
characteristic time L/v t . We assume that e = p/L is a small parameter, with p = v t /fl 
and Q = ZeB /Mc the characteristic gyroradius and the characteristic gyrofrequency of 
the species of interest, and Bq ~ |V x A| the characteristic magnitude of the magnetic 
field. We assume that the characteristic time and length scales in the electrostatic 
potential are the sound gyroradius p s = c s /Qi and the sound characteristic time L/c s , 
i.e., ip(r/p s , c s t/L), where ip has derivatives with respect to its arguments of order unity. 
Here c s = ^/T e0 /mi is the sound speed, rrii and f2j = eBo/rriiC are the mass and 
the gyrofrequency of the dominant ion species, usually singly charged, and T e0 is the 
characteristic electron temperature. The assumption on the scales of the electrostatic 
potential can be easily relaxed to account for other time and spatial scales. Since the 
electrostatic potential ip is a quantity that enters the equations of the different species, 
it is normalized using parameters that do not depend on the species, in particular the 
characteristic electron temperature T e0 , the magnitude of the electron charge e and the 
mass of the dominant ion species m;. This normalization will be useful in Poisson's 
equation, where several species appear. The new, non-dimensionalized variables are 

**=-^,**=.£, (6) 
L L v t B L e s T e0 Mvf 



giving 



with 



£ x (r,v,r,v,t) 



-A(r) + v 



W), (7) 



//*(f,v,i) = i* 2 + Ae0(f/Xe,i/T). (8) 
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Here, T = Mu ( 2 is the characteristic temperature of the species of interest, e s = p s /L is 
the ratio between the sound gyroradius and the characteristic scale length, 



A = — = (9) 
is the ratio between the sound gyroradius and the gyroradius of the species of interest, 



r= V -i = J^L (10) 
c s V T e0 M 1 ; 

is the ratio between the thermal speed of the species of interest and the sound speed, 

and 

yrp / rp \ 3/2 




To \T Q 

Even though the electrostatic potential is small, its perpendicular gradient is not. 
This assumption has been formally implemented by writing h.e(p(i/Xe,t/r). Here and 
in what follows we assume A ~ A ~ r ~ 1. This is the maximal ordering that 
contains in it several interesting regimes as subsidiary expansions in A, A and r. The 
form Ae<p(r /Ae, i/r) is, however, somewhat deceiving because the gradients along the 
magnetic field lines must be small, that is, b • V(p ~ 1 <C 1/Ae, with b(r) := B/£> the 
unit vector parallel to the magnetic field. It is possible to formalize this condition by 
writing the functions in flux coordinates s(f), if>(i) and a(f) such that b = dr/ds and 
B = Va x V?/>. In these variables, the potential is given by 

(p = 0( s (r), V(f)/Ae, a(r)/Ae, i/r). (12) 

To simplify the notation, we will often use (p(i±/\e, fii,t/r) instead of the most complete 
expression in f[T2j) . Where no confusion is possible, we will write <p(r,f). In any case, 
we always assume 

b(r)-V^(r,t)~l (13) 

and 

VrM^i) ■= b(f) x (V f ^(r,t) x b(r)) ~ — . (14) 

At 

Note that e is species-dependent whereas e s = Ae is not. In Sections E] and S] 
where we compute the gyrokinetic phase-space Lagrangian and the equations of motion 
of a single species, e is the natural expansion parameter. However, in Sections [5] and 
El devoted to the gyrokinetic Poisson's equation, we need to consider several different 
species and e s is the appropriate, species-independent small parameter. 

Finally, a notational remark is in order. In Sections E] and H] we will be very careful 
to exhibit the dependence of our results on the mass, charge and temperature of the 
species through the parameters A, A and r. Although at some places this may seem 
unnecessary and awkward (and it would be if our objective were to treat always a single 
species), it is very convenient to write Poisson's equation and the gyrokinetic phase- 
space Lagrangian for a mixture of species in Sections and 

From now on we will drop hats ~ in the normalized expressions. 
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3. Phase-space Lagrangian perturbation theory 

In this section we follow the general strategy of the applications of Hamiltonian and 
phase-space Lagrangian techniques to magnetized plasmas [TBI 113 EDI Ell E2J E3J EU E5] ■ 
We search order by order in the small parameter e for a change of phase-space variables 
such that only one of the variables has fast time dependence. The gyrophase 9 is 
the fast variable that evolves in the gyrofrequency time scale. The rest of the phase- 
space variables (gyrocenter position R, parallel velocity u and magnetic moment /i) 
evolve with the much slower characteristic time scale L/v t . To achieve this, their time 
derivatives dH/dt, du/dt and dfi/dt will be made independent of the gyrophase to the 
order of interest. Had they depended on 9, they would necessarily show rapid time 
fluctuations on top of the more physically interesting slow time evolution. To make the 
time derivatives dH/dt, du/dt and d\ijdt independent of the gyrophase 9, we search 
for a Lagrangian that does not depend on 9 (it will still depend on its time derivative 
d9/dt). We first review briefly how to perform a change of variables in a phase-space 
Lagrangian in subsection 13.11 As we have already announced, the objective is the 
gyrokinetic Lagrangian to order e 2 in general magnetic geometry. Since the calculation 
is quite long and complicated, we have sketched the derivation in subsection !3.2l to offer 
the reader a global perspective of the formalism. This subsection also contains our own 
proof that the algorithm to find the gyrokinetic change of variables can be carried out 
to any order and that there exists an adiabatic invariant /i to arbitrary order. Finally, 
in subsections 13.31 and 13.41 we address the calculation of the gyrokinetic Lagrangian to 



second order in detail. Some of the algebra is relegated to Appendix B and Appendix D 
The results to first order are compared with the iterative method in [17J in Appendix C 
Before proceeding, we must mention that the phase-space Lagrangian (or 
Hamiltonian) approach to gyrokinetic theory has been geometrized [221 E2] • We have 
chosen not to use the language of differential geometry to make the paper accessible to 
a broader audience, without losing mathematical rigor. The reader familiar with the 
geometrical tools will realize that every step of our presentation can be translated into 
that language in an obvious way. 



3.1. Transforming to new phase-space variables 

Consider a transformation T that can be time dependent to a new set of gyrokinetic 
phase-space coordinates {Z a } b a=1 = Z. We write[§| X(Z, t) = (r(Z, t), v(Z, £)) = T{Z,t). 
The phase-space Lagrangian ([7]) can be easily written in the new set of variables by 
using the chain rule, giving 

£ z (z,z,t) = }Z T ^ t ^ - # z ( z >*)> ( 15 ) 

a=l 

§ Note that in part of the literature [25] T stands for the inverse of the transformation that we call T. 
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~A(r(Z,*))+v(Z,i) 



where 

r a (z,t) 

and 

iJ z (Z,t) = tf x (r(Z,t),v(Z,£),t)- 



dr(Z,t) 



-A(r(Z,t))+v(Z,t) 



dr(Z,t) 
dt 



(16) 



(17) 



By finding the stationary points of the action a z [Z(t)] = J^ 1 C z (Z(t),Z(t),t)dt with 
respect to variations of Z(t) subject to the conditions Z(t = t ) = Z and Z(t = ti) = Zi, 
we obtain the new equations of motion 



d fd£ z 



dt \dZ a 



dC 2 
~dZ~ c 



a 



1,2 6. 



Note that the specific form in (j!5p implies that the equations of motion can be written 

as 



dZ? dH z dT r 



a/3- 



dt 



+ 



dt 



a 



1,2,. ..,6, 



with L a/ 3 the 6x6 antisymmetric matrix 



(19) 



(20) 



dZ a dZi 3 ' 

Although our gyrokinetic change of variables has an explicit time dependence due to the 
contribution of the electrostatic potential, we will show that it is possible to choose the 
functions T a such that dY a /dt = 0. In our derivation we impose then that dT a /dt = 0. 
Consequently, we drop the last term in (TT9l and write the equations of motion as 



2> 

13=1 



dZ? 8H Z 



a/3- 



dt 



dZ a 



a 



1,2, 



6. 



(21) 



From expression f[2"Tj) we define the Poisson bracket 



{F,G}= pa 



a,/3=l 



, dF dG 
dZ a dZ? ' 



(22) 



with P Q/3 = (L 1 ) Q/3 the inversq]]] of the antisymmetric matrix defined in (120]) . Then 

dZ a 



dt 



{Z a ,H*}, 



a 



1,2 6. 



(23) 



Noting that L a p satisfies (lA.ip with n — 3, the proof in Appendix A guarantees that 
f[22|) actually defines a Poisson bracket, i.e., that for any three functions Fx, F2 and F$, 
the bracket satisfies skew-symmetry 

{F U F 2 } = -{F 2 ,F 1 }, (24) 

1 1 In Section [4] and |Appcndix E| we show that the matrix L a p of our particular problem is indeed 
invertible. 
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the Leibniz rule 

{F U F 2 F 3 } = {F 1 ,F 2 }F 3 + {F 1; F 3 }F 2 , (25) 

and the Jacobi identity 

{F 2 , F 3 }} + {F 3 , {F h F 2 }} + {F 2 , {F 3 , F 1 }} = 0. (26) 

In general, it is impossible to find a change of phase-space variables that makes the 
Lagrangian as written in (fT5l) independent of gyrophase. However, the time derivatives 
dZ a /dt may be gyrophase independent even if the phase-space Lagrangian is not. This 
apparent discrepancy is easily solved considering that the equations of motion remain 
the same if instead of the Lagrangian (fl5il we employ 

dS 



£(Z,Z,t) = £ z (Z,Z,t) + 



dt 



(27) 



where the function S(Z, t) depends on the phase-space variables Z and t. Indeed, finding 
the stationary points of the action 



a[Z(t)} 



C(Z,Z,t)dt 



to 



£ z (Z,Z,t)dt + 5(Z 1 ,t 1 )-5(Z ,t ) 



(28) 



to 



with respect to variations of Z(t) subject to the conditions Z(t = to) = Zq and 
Z(t = t\) = Z\ gives the same equations of motion as finding the stationary points 
of the action o~ z [Z(t)] because both actions differ only by terms that are held constant. 
Since the Lagrangian is not unique, we are not going to search for new phase-space 
variables Z such that the phase-space Lagrangian C z in (j!5p is gyrophase independent, 
but such that there exists a function S for which the phase-space Lagrangian £ in fT27|) is 
gyrophase independent. This is equivalent to requiring that the time derivatives dZ a /dt 
be gyrophase independent. Thus, we are searching for both the change of variables 
Z and the function S such that the Lagrangian C in f )27|) is gyrophase independent. 
Explicitly, 

6 _ J7« 

£(Z, Z, t) = r«-r- - H(Z, t), (29) 



a=l 



where 

and 
H(Z,t) 



r a (z,t) 



-A(r(Z,t))+v(Z,t) 



dr(Z,t) dS(Z,t) 



tf x (r(Z,t),v(Z,t),t) 



-A(r(Z,t)) + v(Z,t) 



dr(Z,t) dS(Z,t) 



dt 



dt 



(30) 



(31) 



are gyrophase independent. In what follows, Z only refers to the gyrokinetic phase-space 
coordinates. Notice that equations ffT9|) . f )20|) . f l2T|) . f T22l and fT23|) are valid for the new 
Lagrangian C. Simply replace r o by T a and H z by H. 
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3.2. Obtaining the new gyrokinetic variables 

As advanced in the Introduction, the detailed computation to obtain the final form 
of the gyrokinetic Lagrangian to order e 2 is rather involved. This is why we devote 
this subsection to schematically show the steps leading to the determination of the 
change of variables and the function S order by order. We also give an easy proof 
that the algorithm can be carried out up to arbitrary order (although in practice the 
computations would become prohibitively difficult). The detailed calculation to second 
order is done in subsections 13.31 and 13.41 

Our transformation to new phase-space coordinates will be denoted by T e 

(r, v) = T e (R, u, /I, 9, t) = T NP , e Tp, e (R, u, /x, 9, t), (32) 

where R is the gyrocenter position, u is the gyrocenter parallel velocity, /i is the 
magnetic moment and 9 is the gyrophase. For convenience, we have written the 
transformation T e as the composition of two other transformations that we call non- 
perturbative transformation T^p^ and perturbative transformation Tp j£ . 
First, we perform a non-perturbative change of coordinates 



r, v 



7jvp, e (Z s ) = T NPje (R g , v\\ g , fj, g , 9 g ), (33) 



where R g , v\\ g , \i g and 9 g are lowest order approximations to the gyrocenter position, 
parallel velocity, magnetic moment and gyrophase. Their detailed definitions can 
be found in subsection 13.31 For simplicity, we will sometimes use the notation 
{Zg}a=i = Zg = {R g , v fig, 9 g }. Physically, the non-perturbative change of variables 
in subsection 13.31 amounts to saying that to lowest order the gyromotion is circular. This 
is a consequence of the smallness of the gyroradius and the fact that the corrections due 
to the electrostatic potential are of next order. After the change of coordinates we add 
the total time derivative of a function S^piRg, Hg,9 g ) to the Lagrangian. The details 
of the calculation are given in subsection 13.31 The final result is 

£ z *= - e A(R g ) + v llg b(R g ) + e r£> + 6 2 r; 

\ al9 g 



dRg 



where 



+ ("A*. + df> + e 2 rf + -.)^f 7) - HV\R 9 ,v llg ,» g ) - eff« (34) 

H^(Rg,V llg ^g) = ^1+^13^) (35) 



and the rest of the terms are defined in subsection 13.31 Notice that with the notation 
in (I34p we have made explicit the fact that d9 g /alt ~ e -1 , i.e., its time variation 
is of the order of the gyrofrequency time scale and hence much faster than the 
evolution of the rest of the phase-space variables. The change of variables is non- 
perturbative and will give contributions to all orders in e. All the higher order terms 

r R (R-S^HS'/^A)) r R (R-giVWgiVgiOg), •••> r ^ (R g , V\\ g , fJ, g ,9 g ), Vg ) (R g , V\\ g , fl g , 9 g ), ... 

We write t explicitly in (r, v) = T e (R, u, fx, 6, t) because the transformation is in general time 
dependent. 
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and H^(Hg±/e, R g u, fj, g ,9g,t) depend on gyrophase. In the Lagrangian flMj) . the 
Hamiltonian has only the first order correction H^\R g ±/e, R g \\, n g ,9 g ,t), with the 
higher order corrections being exactly zero. This correction H^(R g ±/e, R g \\, fi g , 9 g , t) is 
the only term in the Lagrangian that has strong perpendicular gradients because it is 
the only contribution that depends on the electrostatic potential. 

The gyrophase dependence in the Lagrangian (j34p must be eliminated with the 
definition of the gyrokinetic variables order by order. Since we only calculate the 
gyrokinetic variables to some order, we truncate the expansion in e to the order of 
interest. In general, we need to keep 



1 n 

-A(R 5 ) + v\\ g b(R g ) + cTg (R 9 , v llg , fig, 9 g 



e 

i=i 



dR g 
~df 



+ 



-e^g + J2e i+1 rf(Rg,v ]]g ,ii g ,9 g ) 



i=i 



d9^ 
~~dt 



-H<® (R g , v l]g , fj, g ) - etf « (R g Je, R g]] ,fi g , 9 g , t) + 0(e n+1 , e n+2 ). (36) 

Here, we have not written explicitly d9 g /d(t/e). Instead, we keep the terms that are 
multiplying d9 g /dt to higher order in e. The notation 0(e n+1 , e n+2 ) will be extensively 
used in this paper and indicates that the terms of order e n+1 that we have neglected are 
either proportional to dR g / dt or are in the Hamiltonian, and that the terms of order e n+2 
that we have neglected are proportional to dv\\ g /dt, dfi g jdt and d9 g /dt (the Lagrangian 
fl36]) does not contain terms proportional to dv\\ g /dt or d/j, g /dt). When we perform the 
expansion to obtain the gyrokinetic variables order by order, it will be apparent that 
this notation is convenient because we need to keep some terms to 0(e n ) and the rest 
to 0(e n+1 ). 

Employing expression ( I36p . we find the gyrokinetic variables by eliminating the 
gyrophase dependence order by order. For now, it is enough to simply write the 
transformation 

(R g , v\\ g , fig, 9 g) = T P>e (Z, t) = T P>e (R, u, n, 9, t) (37) 
to the order of interest as 

n n n n 

R g = R + Y^ e m R m , v\\ g = u + ei u>i, N = fi + 6 9 = d + J2 ( 38 ) 

i=l i=l i=l i=l 

In subsection 13.41 we show the connection of this expansion with Lie transforms. 
Sometimes we will use the abbreviated notation {Z a }^ =1 = Z = {R, u,fi,9}. Notice 
that the variable R is found to an order higher than the rest. The corrections R n +i, 
u n , jl n and 9 n are obtained by imposing that the Lagrangian is gyrophase independent 
up to terms of order 0(e n , e™ +1 ), where the terms of order e n are in the Hamiltonian or 
are terms proportional to dR/dt, and the terms of order e n+1 are terms proportional to 
du/dt, dfi/dt and d9/dt. To demonstrate the procedure, we show schematically how to 
obtain the first corrections R2, u±, fi\ and 9\. We then argue that the same formalism 
can be extended to arbitrary order. The proof presented here demonstrates that the 
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gyrokinetic variables can be consistently calculated order by order without running into 
problems. Recently, this has been a controversial issue [3"9~ ] I4TJ ] |4"T]. 

To calculate R 2 , %, jli and 9\, wejieed to express the Lagrangian ( 1361) in the new 
gyrokinetic variables to 0(e, e 2 ), givin 



£ 2 



-A(R) + ub(R) + e ( R 2 • VrA(R) 



Mib(R) + r{J ) (R,w,/i 



+ 



eA(R) + t 2 uh{R 



dR 2 
~dt 



e/i + e z ( — pL\ + Tg 1 ^ (R, u, n 



d9 
~dt 



dR 

~dt 

- e 2 /i- 



d6x 
~dt 



//i5(R) + (R ± /e, % (i, 9, t)] + 0(e 2 , e 



(39) 



Notice that in the functions r^ } (R 9 , v \\ g , n g , 6 g ), T d L) (R gj v\\ g , fi g , 6 g ), if (0) (R 9 , v\\ g , /i g ) 
and H^\R g ±/e, R g \\, n g ,9 g ,t) the variables R g , Vu g , [i g and 9 g are replaced by the 
gyrokinetic variables R, u, /i and 9. For example, from the definition of in (|35|) . 
we find 

# (0) (R, u, //) = )V + /i£(R) . (40) 

The correction uu\ + /2iS(R) + H^(R±/e, R\\, fi, 9, t) to the Hamiltonian contains both 
the correction and the result of Taylor expanding H^°'(R g , v\\ g , n g ) = H^°\R + 
e 2 R 2 + . . . , u + eu\ + ...,/! + efx\ + . . .) around R, u and /i, i.e., 



Ur- 



du 



+ jl 1 ^— + H (1 \R ± /e,R ll ,fx,9,t) 



dfx 



uu 



x + hB{R) + (R±/e, R\,fi, 9, t). 



(41) 



In subsection 13.41 we show that the expansion around R gives a term that is an order 
higher and can be ignored to this order. 

As indicated in subsection 13.1] we can always add a time derivative to the 
Lagrangian f l3"9~|) . By doing so, we will get a form of the Lagrangian in which the 

(2) 

derivatives of a function S P (R±/e, R, u, /i, 9, t) enter. By imposing that the Lagrangian 
is of a specific form, we first obtain the corrections R 2 , Ui, pL\ and 9\ as functions of 



Sp^ and we then find Sp' by integrating a simple differential equation. We add to the 



7(2) 



Lagrangian (1391) the total time derivative 
d_ 

dt 



e 2 S ( p ] (R±/e, R, u, n, 9, t) - eA(R) • R 2 - e 2 wb(R) ■ R 2 + e 2 /A 



(42) 



(2) 

where we have taken into account that S p depends on R in two different ways: a fast 
dependence due to the potential that has wavelengths on the order of the gyroradius, 
and a slow dependence due to the background magnetic field. Notice that adding the 
time derivative (142]) to the Lagrangian ( 139]) eliminates all the terms proportional to the 

+ Observe that in principle, the term e 2 A(R) • dR^/dt should be included in C z (|55]l . However, adding 
the time derivative of — e 2 A(R) • R3 eliminates any R3 dependence to this order. The same can be said 
about (J33J) and (|114j) . where the terms e" +1 A(R) • dR n+2 /di and e 3 A(R) • dR 4 /eft could be added. It 
is very easy to show that the algorithm gives, of course, exactly the same results with either choice. 
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C 



^A + M b + e(BxR 2 + M 1 b + rJ ) + V (R±A) ^ 



(2) 



dR 

~dt 



du 



du J dt 



+ e 2 \e 1 + 



OS 



(2)~ 



dp, 
~dt 



-en 



r 



(1) 



dS 



(2)~ 



06 



d9 
~db 



H (o) 



-e(mi 1 +/ii5 + # (1) ) +0(e 2 ,e 3 ), (43) 
where we have used 

R 2 • V R A - V R A • R 2 = (V R x A) x R 2 = B x R 2 , (44) 
trivially deduced from the identity 

r)A- r)A- 

[V R A - (VrA) t ]^ = - — 1 = J> ife (V R x A),. (45) 

* i k=i 

Here is the Levi-Civita symbol and the superscript T stands for matrix transposition. 
Notice that in equation (|43p we have not made explicit the dependence of the functions 
on the gyrokinetic variables, but it is assumed that all the terms are functions of R, u, 
\x and 6. By imposing that the Lagrangian (T4"3"|) be equal to 



-A + uh + eT l £ 

e 



dR dO 



dt 



dt 



where T^, 



H 



H( \R, u, //) and are gyrophase independent, we obtain the 



(0) 



f(i) 



eW ±J + 0(e 2 ,e 3 ^ 



(46) 



equations 



R 2 



dS 



(2) 



du 



-b(R) + 



B(R) 



b(R) x 



rg ) (R ) tt >A i,e)-r^ + v C a ±/e)1 s^ 



:(2) 



"1 



b(R) 



r^-rg } (R, w 



aa(2) 

/ i 1 =ri 1) (R,n, /U ,^) + ^ 



and 



OS 



(2) 



(47) 
(48) 

(49) 
(50) 



(2) 

The corrections R 2 , Ui, jli and #1 can then be found if we obtain S p . To do so, we use 
that 

uu 1 + B{R)fi 1 + H {t \R L /t,R\ b ^9,t) = H {1) . (51) 

First, we take the gyroaverage of this equation to obtain H 1 . Using the results in 
equations and fT4"9l . the gyroaverage of equation flBTj) becomes 



(i) 



ub(R) 



R 



r R } (R, M ,/i, 



s(R)(rW(R,«,/i 



+ (H^(Rje,R ll , f i,e,t)) 



(52) 
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where (...) = (2n) 1 § d6 (. . .) is the gyroaverage holding R, u, fi and t fixed. Notice 



K 1 ) 



that we have been able to obtain if vv without explicitly finding Sp\ Solving for Sp 



7 (2) 



is easy once is known. Combining equations ( 151]) and (|52|) . we find 



95 



(2) 



■it 



89 



B(R) 



-b(R 



S(R) 



(R, u,fi,9) — (r$ (R, u, 
(R, u,fi,9) - (rf) (R, m, < 

1 [HW(R ± /e,R\\,ii„6,t) - (H^(R ± /e,R\\,fj,,e,t))] 



(53) 



It is then straightforward to obtain S p by integrating in the gyrophase 6. The specific 

T7« 



calculation of the corrections R2, Hi, #1 and H is done in subsection 13.41 Here 
we only want to demonstrate the procedure. Notice that the final solution depends on 
our choice of - and in general on our choice of the final expression for C. However, 
once the choice is made, is completely determined. 

Once R 2 , ui, ftt, 8~i, and Sp^ have been determined, we can prove by induction 
that the procedure gives the corrections to any order. Assume that the corrections R»+i, 
Ui, fli, 8~i and H % and the functions Sp +1 ^ have been obtained up to % = n — 1 in such 
a way that the Lagrangian is of the form 



n— 1 

-A(R) + ub(R) + J^rf 



dR 

~dt 



e/i 



d6 
dt 



H 



(0) 



n-l 



i=l 



+0(e n ,e' 



n n+l N 



(54) 



Then it is possible to obtain the corrections R n +i, u n , jl n and 9 n , the function 4 +1) and 
the phase-space Lagrangian to order 0(e n ,e n+1 ). To do so, the Lagrangian is written 
to 0(e n ,e n+1 ). In general, we cannot guess the exact form of the Lagrangian to such 
high order without doing the calculation order by order, but the terms that contain the 
corrections R„+i, u n , fl n and 6 n are very easy to obtain. We find that 

n— 1 



c 2 



1=1 

1 

€ 



+ 



+ 



J t (e 4+1 4 +1) - e l A ■ R i+1 - e i+1 ub ■ R l+1 + e i+l ^9 

n-l 

A(R) + ub(R) + + e ™ ( R «+i • VrA(R) + uJ>(R) + f £ 5 

e n A(R) + e" +1 wb(R) 



8=1 



dR 

~dt 



-e/i + e 



n+l 



-u + r y 



(n) 



1 n+lf(n)_ 1 n+lWn) ""P 
U dt " dt 



r n+l d®n 
^ dt 



H 



(0) 



n-l 



(i) 



i=l 



+ /x n S(R) + 



0(e n+l n+2N 



(55) 



Here we have just separated the terms of order 0(e n ,e n+1 ) into those that depend on 
the corrections R„+i, u n , jl n and 9 n , and the rest that we have lumped into the terms 



f,M f,(n) p(n) p(n) 



R 5 



J- u 



and #H 
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The form of the Lagrangian f l55|) is very similar to the Lagrangian f l39|) . We can 
then use the same procedure. We add the time derivative 

d_ j e n+l£(n+l) ( Rj y gj R ^ ^ ^ f) _ £ „ A(R) . R ^ _ ^+1^^ . + gH+l^j (gg) 

to cancel all the terms that are proportional to the time derivatives of the corrections 
R n+ i, u n , Ji n and 6 n . Then, by imposing that the Lagrangian ( 155]) plus this time 
derivative be equal to 



Z = ( i A(R) + «*R) + £ ^ ) ■ ^ - *i - 77'"' - £ e -H 



(57) 




BfR' 



■b(R)x(fS ) -rS ) + V (R±A) 4 n+1) )> (58) 



« n = b(R)-(rS ) -fL n) ), (59) 



and 



An = f W + (60) 



flp(n+!) 

ft, = -f<"> - ™E—. (61) 



The n-th correction to the Hamiltonian becomes 



H (n) = „b(R) • (r<? - (f W)) + 5(R) (f W } + (#<»>) , (62) 



and the equation for Sp 1-1 " 1 ^ is 



#(»)-/#(*))). (63) 



B(R) 

In subsections 13.31 and 13.41 we obtain the phase-space Lagrangian to 0(e 2 ,e 3 ) 



Specifically, in subsection 13. 31 and Appendix B we derive equation (134"]) . In subsection 13.4 



and Appendix D we use the perturbation procedure explained here to go from equation 



341) to the final result. 



3.3. Non-perturbative change of variables 

We perform a change of variables (r, v) = T NP)€ (Z g ) = Tjvp, 6 (R g , v \\ g , fi g , 9 g ) defined by 
r = R g + ep(B. g ,fig,9g), (64) 

and 

v = v llg b(R g ) + p(R g , Hg,9g) x B (R fl ) , (65) 
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with the gyroradius vector defined as 

p(Rg,fXg,9g) = " 



2^ 



BR n 



[sm9 g ei(R g ) - cos6» 5 e 2 (R 9 )] . 



(66) 



The unit vectors ei(r) and e2(r) are orthogonal to each other and to b = B/JB, and 
satisfy e x x e 2 = b at every location r. Physically, R g is the guiding center position, v\\ g 
the velocity parallel to the magnetic field at the guiding center position, \i g the lowest 
order magnetic moment, and 9 g the lowest order gyrophase. For a homogeneous static 
magnetic field and in the absence of electric field, the change of coordinates defined 
by (|65|) and ( 166]) exactly eliminates the gyrophase dependence. Note in passing 

that it is a well defined change of coordinates. These formulae explicitly give {r,v} 



as a function of Z 9 = {R s , v\\ g , fi g , 9 g } and it is easy to see that the transformation 



is invertible for small e: it is clearly invertible for e 
continuous in e. 

We substitute the relations X(Z 



and the transformation is 



g j, given in (1611) and (1551) . and 



dX a 
dt 



Zg, Zg) 



j2 dxa ^ *? 



azS 



9 ■ 



a 



1,2, 



(67) 



into the non-dimensionalized Lagrangian (0), £ x (X,X,t). The resulting La grangian 
that we denote as £ X (X(Z 9 ), X(Z 9 , Z g ),t) differs from the Lagrangian £ z « in (JMD and 



(185|) by the time derivative of a function Snp and even though both Lagrangians give 
the same equations of motion, we have decided to stress the difference. The Lagrangian 



£ x (X(Z 9 ),X(Z 9 ,Z 9 ),t) is 



(X(Zg),X(Zg,Zg),t) 



~-A(R g + ep) + v\\ g b g + p x B 5 



d_ 

dt 



[Rg + ep) 



-HW-eHM, 

with H^°'(R g , v\\ g , n g ) defined in (155|) and 

i/ (1) (R 9± /e, ^||, ^, t) = A(0)(R 9± /Ae, R glh fi g /X 2 , t/r) 
+A0(R 9± /Ae, R glbf i g /\ 2 , 9 g + tt0(-A), t/r), 



(6* 



(69) 



where 0(x) is the Heaviside step function, with @(x) = 1 for x > and 0(x) = for 
x < 0. Any magnetic quantity with subindex g is evaluated at R g , e.g., B 5 := B(R S ). 
We write p = p(R 9 , when no confusion is possible. 



We have defined a new function 
between and tp) given by 



[R g , fj, g ,9 g ,t) (notice the difference in the font 



(Rg, fig, 9g, t) \= (f(Rg + €p(Rg, fig, 9 ' g) , t) . 



Then (0) is the gyroaverage of 



((j))(Rg,Hg,t) 



1 

2~ix 



2tt 



d9g 4>(Rg, flg,9g,t), 



(70) 



(71) 



Phase-space Lagrangian derivation of electrostatic gyrokinetics 



18 



and 4> the gyrophase dependent piece, 

(fr(Rg,llg,9g,t) = <j)(Rg, fig, 9g, t) ~ ((/)) (Rg , fig , t) . (72) 

We now prove that the notation in fl69|) . 

0(R 9 , fig, 9 g , t) = 0(R 9± /Ae, R gll ,fi g /X 2 , 6 g + 7r6(-A), t/r), (73) 

is appropriate. First, we show that fi g is always divided by A 2 and that the sign of A 
determines the phase of 6 g , and later we demonstrate that if conditions ffTBl and ffT4l) 
are satisfied, then 



b(R s )-V Ra <XR fl ,/i fl ,M)~l 



(74) 



and 



V R9X 0(R 9 ,/i 9 ,^,t) := b(R 9 ) x (Vn a <l>(R B ,n g ,e g ,t) x b(R 9 )) ~ — (75) 

are also satisfied. To show that fi g always appears divided by A 2 and that we need to 
add 7r to 6 g when A is negative, it is enough to realize that depends on ft g and 6 g 
through the dependence of ip on r/Ae = R/Ae + p(R 9 , fi g ,9 g )/ A and that p as defined 
in (16"6"j) only depends on fi g through the multiplying term y /TI g '- It is then obvious that 
r/Ae = R g /\e + p(R g ,fi g /\ 2 ,6g + nQ(-X)). To prove that (Q3]) and (JUD imply (J71J) and 
( 175]) . we employ 

Vr>(R 9 , = V Rg v3(R 9 + ep(R g , ff ), t) 

= V r <^(r, t) + eV Kg p(R g , fi g , 9 g ) ■ V r y?(r, t), (76) 

with cVr 9 p ■ V r (/? ~ cVr 9 p ■ V r± y? ~ A -1 ~ 1. Then, using equations ( 1X3]) and (IH]), 
it is easy to see that equations ( 1T1"|) and f lT5"j) are correct. Note that when we write 
H^(R g ±/e, R g \\, fi g ,9 g ,t) in (169]) we are emphasizing the dependence on e because the 
asymptotic procedure is based on expanding in e -C 1. The dependence on A, A and r 
is only written explicitly in the function <p. 

We now show how to simplify f )68|) . Employing 



2B„ 



and 



Vr 9 P 



dp J_ 

2flg 



dp , 



P - (V R9 b 9 • p)b 9 + V Rg e 2 3 • e lg [p x b 



9)1 



(77) 
(78) 

(79) 



we write the Lagrangian in 



as 



£ X (X(Z,),X(Z 5 ,Z 9 



),t) 



-A(R g + ep) + U|| fl b g + p x B g + V Rg p ■ A(R 9 + ep) + 



e( 2/J s V R9 e 2s ■ e lg - W|| 9 V R9 b 3 • p 



rfR 



+ — A(R 3 + ep) ■ p 



dt 



A(R g + ep) 



<9p 



2e/i g 



(80) 
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To obtain (JZZJ) we have used V Rg ei 3 = -(V Rg b ■ e lg )b g - (Vr^j ■ e lg )e 2g and 
V Rg e 29 = -(V Rg b • e 2g )b g + (V R9 e 2g ■ e lg )e lg . To simplify the Lagrangian flgO]), we add 
the time derivative of 

dfj,' 



SNp(Rg, /lg, 9g 



2/i 



f A(R 9 + ep(R 9 , n' g , 6 g )) ■ p(R g , ^ 6 g ). (81) 



As a result we find 



-A(R 9 + ep) + v\\ g h g + pxB g + V Rg p ■ A(R 5 + ep) + V Kg S NP 



+ 



e^2/i 9 V Rg e 29 ■ ei s - ^|| 9 V Rg b s • p 



c/R n 



dfl g 
~dt~ 



12) 



In Appendix B we prove that 
1 1 

V Rg SV P = — -A(R g + ep) + -A fl - p x B g - V Rg p ■ A(R g + ep) - e/i 5 V Rg e 2g • e lg 
\p' x [B(R g + ep') - B g ] + e[(p' x b g ) ■ B(R g + ep')]V Rg b, ■ p' 









Jo 




By 


and 






dp 




de g 



2/i 



+ -p[b 9 • B(R 9 + ep') - £,]V Rg e 2 , • e lg 



A(R 9 + ep) + e/2 g + — / d/i g [b s • B(R 9 + ep') - £ 



(83) 



34) 



9 JO 



where we use the abbreviated notation p' = p(R 3 , [i! gl 9 g ). Substituting equations (J83 
and (l84"j) into the Lagrangian (1821) finally gives 
1 



with 



-A 5 + D|| ff b g + eAT R 



^ + (-e N + e 2 Ar e ) ^ - H® - effW, (85) 



AT R = /i 9 V Rg e 29 ■ e lg - W|| 9 V Rg b 9 • p g 



"» da' { 1 

£f -p'x[B(R g + ep')-B,] 

o ¥0 I e 



2u' 

+ [(p' x b 9 ) • B(R, + ep')]V Rg b, . p ' + J^[b g . B(R g + ep') - BJVh,^ • e lg j>(86) 



and 



Ar 6 



1 f^ 3 



eB 



dfj,' [b g ■ B(R 9 + ep') - B, 



37) 



5 ^0 



It is easy to write the Lagrangian (1851) order by order. We use 

J2 



B(R g + ep) = B, + ep • V Rg B 9 + -pp : V Rg V Rg B, + 0(e 6 ), (88) 



where our double-dot convention is ab :M= b- M -a, to obtain 



R 



+ [-e, g + e^ + e^) d A 



(89) 
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where 

r R = ^9 V R 9 e 29 • e lg - v || 9 V Rfl b 9 • p - -(p • V R9J B 3 )p x b 9 



+^-VR 9 b 9 -(pxb s )]B 9l (90) 

1 B 
r R = qPP '■ Vr 9 Vr 9 B 9 X p - -^[p • V Rg b 9 • (p x b 3 )]V Rff b 9 • p 



2^_ 

35 s 



p-V Rs B s )VR 9 e 29 -e l9l (91) 



r (l) _ 


2Hg 




3B g 


and 




F (2) _ 






45, 



pVn Q B g (92) 



PP . V R9 V R9 B 9 b r (93) 
5.^. Perturbative change of variables 

In this subsection we find a new set of coordinates {R, u, fi, 9} that makes the Lagrangian 
(I59"|) gyrophase independent. We employ the procedure described in subsection 13. 2[ The 
transformation (R,, v\\ g , 9 g ) = Tp i£ (R, u, /i, 9, t) is customarily written in the form of 
a Lie transform 



T Pj£ = T X T 2 T 3 . . . , (94) 

where 

T n = exp e n+1 R n+1 -Vn + e n (u n —+fi n — + 9 n —j . (95) 

Instead of this form, we use the form in (13"5|) that we find more convenient. The 
connection between the two arrangements is trivial. To first order we find 

R 2 = R 2 , Ui = «i, fi\ = 9\ = 1 . (96) 

To second order, the relation is 

1 ui <9R 2 /i! <9R 2 9i <9R 2 

-R 2 ■ V (RxA) R 2 + y — + y — + - — 



R 3 = R 3 + -R 2 ■ V (RxA) R 2 + -^-^ + + (97) 



and 



, 1 r» V7 uxdui (Midu! 9 1 du 1 

u 2 = u 2 + -R 2 ■ V(R ±/e) «x + y ^ + y ^ + Yl"' ( 98 ) 

^ = ^ + 2 R2 ' V ^/^ + Y ^7 + T^7 + 2 W (99) 

0, = « 2 + ^ • V (Ri/<)9l + ^ + f ^ + (100) 

In subsection 13.4.11 the corrections R 2 , u±, /ii and 9i are calculated following the 
procedure in subsection 13.21 and the Lagrangian is obtained to 0(e,e 2 ). In subsection 
13.4.21 the Lagrangian is obtained to next order. It is possible to do so without explicitly 
obtaining R 3 , u 2 , fi 2 and 9 2 . 
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3.4-1- Perturbative change of variables to first order. We obtain the first-order 

T7« 



gyrokinetic correction to the Hamiltonian, H , by employing equation (|52|) . We need to 
know that (if^R, u, /i, 0)) = 0, (#«(R ± /e, % //, 0, r)) = A(0)(R ± /Ae, % /i/A 2 , t/r) 
and 

(r£> (R, it, /i, 0)) = /xV R e 2 -ej + ^bx V R 5 - ^hh • V R x b, (101) 



where we have used that 
<PP> = f (I 



bbl 



(102) 



with I the unit matrix. For the remainder of the section, whenever we do not write 
explicitly the arguments of the functions, it will be understood that they are evaluated at 
(R, u, fi, 8), i.e. b = b(R), (0) = (0)(Rj_/Ae, R\\, /i/A 2 , t/r), and so on. Substituting the 
values of (r^(R, u, fi, 8)), (r^(R, u, fi, 8)) and (H^(R±/e, R\\, /i, 8, t)) into equation 
(E2D, we find 

~' ~ ""' (103) 



H W = A(0) + uh ■ r£ } - ufjb ■ V R e 2 ■e 1 + ^b-V R xb. 



Notice that we have the freedom to choose as we wish. Our choice will affect the 
corrections R 2 and u±, and the final form of H 1 . To coincide with previous derivations 



in the literature 



we choose 



// 



bb • V R x b 



giving 



H 



(i) 



A(0). 



(104) 



(105) 



In equation (I104p . we have chosen = /iV R e 2 • e x + . . . instead of T-^ = /xbb • Vr^ • 
ei + . . . to manifestly show that the equations of motion are independent of the choice 
of ei and e 2 [4"2] . 

The function S P can be obtained by solving equation fl53|) . Substituting equations 
(EHD, ([90]) and ([92]) into fl53D gives 



# = -B b ' V " b '" + 4 



where we have used that 
PP - (PP) 



1 



P(P x b) + (p x b)p 



PP - (P x b)(p x b) 



: V R b 



2/i 
3~B 



p - V R £? 



A0 

IF' 



(106) 



(107) 



Integrating equation (I106p in the gyrophase gives 



S 



(2) 



U 

~B 



u 



■■ V R b 



b • V R b . (p x b) - - pp-(px b)(p x b 

O - 

2u . , _ _ A$ 

where the function $ is the integral 
$(R ± /Ae, % ///A 2 , 8 + nG(-\),t/r) 

tZ0' 0(R±/Ae, R\\,fM/\ 2 , 8' + 7r6(-A), t/r) 



(108) 



(109) 
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such that ($) = 0. Here we have used that p = d(p xb)/d9 and p(p x b) + (p x b)p = 

-(l/2)d[pp-(pxh)(pxh)]/d9. 

Using S P in the expressions fj47|) , (PfBj) , (|49i) and (1501) , the first order corrections to 

the gyrokinetic variables become 

2u - - ~ 1 - r 

R 2 = --gbb • V R b • (p x b) - -b ^pp - (p x b)(p x b) 



u 



■ V R b - — b x V R b ■ p 



rb x V (R±/Ae )$, 



5 



«i = uh • V R b • p — — p(p x b) + (p x b)p 



: V R b 



2 

/ ' i -— b-V R b p+ 7 p(pxb) + (pxb)p 



and 



: V R b 



Tl A 1/ r 

0i = 7T^ b ' V R b • (P x b ) + F" PP - (P x b )(P x b ) 

8/i L 

1 / fx ^ r, A d ^ 

+ ^ Xb) ' VRB + ^W 



A0 
IT 



: V R b 



(110) 
(111) 

(112) 



;ii3) 



In |Appendix C| we show that this result is equivalent to the result obtained with the 
iterative method in [Tj 



3.4-2. Perturbative change of variables to second order. In this subsection we apply 
the change of variables ( 155|) to the Lagrangian ( 159]) to 0(e 2 , e 3 ). The idea is to write an 
expression similar to ( 155|) with n = 2 so that we can use the technique demonstrated in 
subsection 13.21 

To 0(e 2 ,e 3 ), the Lagrangian f )89|) becomes 



-A + wb + e ( rff + R 2 • V R A + Ul b ) + e 2 I T^ + R 3 • V R A + w 2 b 



^(2) 



<9r (1) 

+ wR 2 - V R b + Ml — ^~ 

ou 



arg> ar£> 



+ 



eA + e 2 ?ib + e 3 ( + R 2 • V R A + Ul b 



R 

dR 

' ~dt 

dR 2 



+ 



-efi + e 2 (-/i! + r 
- e 2 /i + e 3 (-/xi + T 



(i) 
(i) 



-£ 2 + rf+/x 



dt 



dfj, 



+ |^e 2 A + e 3 ub 

9r (i) 
+ 0v u 



dR 3 
~df 



09 



d9 
~dl 



d9 l dfh _ H{0) _ e , + + R{1) s 
dt 



dt 



e 2 (uu 2 + ^+j2 2 B + //Ha • V R 5 + R 2 • V (RxA) # (1) + /ii 



+ 9i 



09 



+ 0(e 3 ,e 4 ), (114) 



where we have used that dT e jdu = 0, that dH^/du = and that T R and T g only 
depend slowly on R. 
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We have seen in subsection I3.2l that to first order we need to add to the Lagrangian 
f)114p the time derivative 
d 



dt 



e 2 S ( p ) - eA R 2 - e 2 ub ■ R 2 + e 2 /A 



(115) 



giving as a result the Lagrangian 
d 
Jt 

-A + ub + eT^ + e 2 (r 3 ■ VrA + u 2 b + f ® 



£ z + ^ ( e 2 ^ 2) - eA R 2 - e 2 uh ■ R 2 + e 2 /A 



^ + (f 2 A + r :! ,/b 



dR, 



, f 3f(2)^ , 3f(2)^M 



e/i + e 3 ( -# 2 + f (2) 
uuz + faB + H® + 0(e 3 ,e 4 ) 



dt dt 



(116) 



with 

f R = r R + m (Vr x b) X R 2 + Ml 



ar 



(1) 

R 



(9u 



ar 



(i) 

R 



a F (l) 

^^ + Vr4 2) 



06 



v (r±/£) r 2 • ( r R } + r 2 • v r a + Ul b) + [ -m + r^ iJ ) v (Rx/e) ^, (117) 



f(2) 



+ R 2 ■ VrA + Wl b 



f (?) = f r( 1} + R 2 • V R A + Ul b 



0R 2 

Ou 

0R 2 



+ (-^ + if ) 
+ f-Aii + r? } 



Ou 
06 1 



(118) 
(119) 



'-.(2) _ r( 2 ) 



rr + A*i- 



dY 



(i) 



9/i 



+ 0i 



ar 



(i) 



a# 



+ r 



R 



R 2 • V R A + uib 



dR 2 

06 



and 



#( 2 ) = ^ + ^R 2 . V R 5 + R 2 ■ V (RxA )# (1) + + #r m 



OH^ „ di/W 95 



(2) 



at 



(120) 



(121) 
(122) 



Notice that in equation ( I117P we are using 

R 2 • V R b - V R b • R 2 = (V R x b) x R 2 , 

(2) 

and we have taken into account that S P depends on R in two different ways. On 
the one hand, there is the dependence on Rj_/e, that was the only dependence that 
was taken into account in subsection 13.4.11 This dependence gives the strong gradient 
V(R ± / e )5'p = — (A/AS)V(R ± /A e )$. On the other hand there is a slow dependence on R 
that gives the gradient 



(2) 



R 



— b • V R b ■ (p x b) - - ( pp - (p x b)(p x b) ) : V R b 



-|(pxb)-V R B 



fv R B - ^bb ■ v R $. 



(123) 
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In subsection 13.21 we showed that by adding the total time derivative 

4- (e 3 S ( p ) (R ± /e, R, u, /I, 9, t) - e 2 A • R 3 - e 3 ub ■ R 3 + e 3 fi9 2 
at V 

to the Lagrangian (11161) and making the result equal to 

z 



24 



(124) 



^A(R) + ub(R) + eT^ 



^1 _ & _ F (0) - eH (1) - e^ 2) + 0{e\ e% (125) 



T(2) 



where we have explicitly set T R = 0, we obtain the equations for the corrections 



Ra 



f (2) , O&P 

u ^ du 



b(R) 



B{R) 



(3) 



b(R)x r R ' + v (R±A) ^ 



? i 2 = -b(R)-fg ) , 



A, = i?> + 



95 



(3) 



00 



and 



02 



p(2) _ aD P 



The correction to the Hamiltonian is 

H {2) = -ub(R) ■ (f®) + B(R) (if \ + (H^\ 



(126) 

(127) 
(128) 

(129) 
(130) 



and the equation for Sp^ is 



OS 



(3) 



09 



" -MR) • ff <? 



BR) 



f (2) 
1 R 



=.(2) 



',(2) 



B(R) 



(3) 

In this article we do not solve for S P and hence we are not able to write explicitly 



Tr(2) 



the corrections R3, u 2 , fi 2 and 9 2 . We only obtain explicitly the correction H . In 



Appendix D we evaluate equation (I130p in detail. The final result is 

H {2) = A 2 ¥*\Rj\e, R, /i/A 2 , t/r, A) + A*g(R ± /Ae, R, w, /i/A 2 , t/r, A) 
+*g ) (R, M ,/i), 



«9<0 2 > 



with 






1 ( 


2A 2 5 2 N 




- — < 


\B V 




1 




4AS 



V(R X /A £ )$ • b X V (R±/ Ae)0 ) - 



7 (R ± /A e )0 x b ■ V R b ■ p 



2X 2 B <9(/x/A 2 
/' 



(132) 
(133) 



2A5 2 



V 



(RL/Ae)< 



pp-( P x b)(p x b) 



V R fi- V {R±/Ae) (0) 

1 



5 



:V R 5 
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U — s 
+ — Vnb : 



p(p x b) + (p x b)p 



4A 2 R ' \d(fi/X- 

P(P x b) + (p x b)p 



and 



(2) 
B 



3?/ 2 /y ~ n 2 <-> 

-^b ■ V R b ■ V R 5 + i^(I -bb) : V R V R B ■ b 
+ ^V R b:V R b + ^ 



(134) 



|V R± 5| 2 



+ 



2B 
3u 2 fi 



U 2 fJ, 

~4B 



v R b : (v R b) T 



u 



2B 2B 2 



b- V R b| 2 + 



u 2 [i fi 2 
~8B ~~ 16 



3/x 2 
45 2 

^ + ^] (V R -b) 2 
85 16 ' 1 R J 



(b • v R x b) ; 



(135) 



where M is the transpose of the matrix M. The final phase-space Lagrangian is given 
then by ( Il25p . We can write it explicitly as 

"1 



-A(R) + ub(R) - e//K(R) 



dH d6 __ „ , 
_- £/1 --/, + ( £ V), 



(136) 



where 

H=^u 2 + (iB(R) + Ae(4>) (RjXe, % /i/A 2 , t/r) + A 2 e 2 ^ 2) (R ± /Ae, R, fi/X 2 , t/r, A) 



+Ae 2 ^ 2 ] (R_l/Ac, R, u, fi, 11/ A 2 , t/r, A) + e 2 ^> (R, u, //) 
K(R) = ib(R)b(R) ■ V R x b(R) - V R e 2 (R) ■ e x (R). 

(2) 



^T,( 2 ), 



and 



(137) 



(138) 



In previous work [24"| [25] . only the contribution was kept because the terms that 
contained the function <fi were assumed to be larger. With the more natural ordering 
(JTJ) , we find the new contributions and W B , demonstrating that magnetic geometry 
and electrostatic potential appear together and cannot be separated. In Sections HJ [5] 
and E] we show that and W B ' modify both the equations of motion and Poisson's 
equation. 

We end this section pointing out that we chose the final form of the Lagrangian 
(HH, with IV = 0, to have the same Poisson brackets as previous authors [25]. There 
are other possible choices, e.g., making the second order correction of the Hamiltonian 
independent of the parallel velocity, condition that can be achieved by defining the 
appropriate T^. 



4. Gyrokinetic equations of motion and Vlasov equation 

The equations of motion are given by (|2"3"|) . To obtain them explicitly we need to find the 



Poisson bracket (1221) that corresponds to the Lagrangian (11361) . Employing Appendix E 
we find the Poisson bracket to be 

fFC\- 1 ( 9FdG dFdG \+ B * ( v * F dG 9F v * r 
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+^F-(bxV^), (139) 
B \\ 

where 

B*(R, u, y) := B(R) + euV R x b(R) - e 2 iiV R x K(R), (140) 
BUR, u, fi) : = B*(R, u, fi) ■ b(R) 



B(R) + eub(R) ■ V R x b(R) - e>b(R) ■ V R x K(R) (141) 



and 



V R :=V R -K(R)^. (142) 
Employing the Poisson bracket in f)139p and the Hamiltonian in f)137p . we find 
• dR ( k 2 d^f B 2 «9# (2) \ B* 1 . f „ n Ae„ 

+ ^V (Rx/A ^r + ^V (Rx/Ae) *g> + A 2 e 3 V R < 2 > + Ae 3 V R *g 

+ e 3 V R ^g ) \ (143) 



u = ^ = -JjB* • V R 5 - Aeb ■ V R (0) - A 2 e 2 b • V R v^ 2) - Ae 2 b • V R ¥g> 

- e 2 b ■ V R ^g } -^[»bx(b. V R b) - eyU (V R x K) ± ] ■ V (R±/A£) (0) 

+ ^V (R±/A ^J) + ^V^*g + A 2 e 3 V R < 2 > + A^V H *g 

+ e 3 V R v&g ) ^, (144) 

and 

a _ d9 1 A 0(0) A 2 e <9^ 2) 03$ Ae ggg gjjff 

dt e X 2 d(/j/X 2 ) X 2 d(/j/X 2 ) € dn X 2 d{n/X 2 ) 6 0// 

tx + Ae 2 — ^ + e 2 -^- - — (K X b) • Uv R 5 



B* V du du J B* } 



+ y V (Rx/Ae) (0) + ^!v (Rx/Ae) *f + ^V (R±/Ae) *g> + A 2 e 3 V R vI> 



+ Ae V R \&0^ + e 3 V R ^ J J . (146) 

Note that we have emphasized the fact that the dependence of the functions (0), 
and \I/^g on R and /x (recall equation ( I137P ) can be fast or slow. For this reason we 
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distinguish between derivatives with respect to the argument Rj_/Ae and derivatives 
with respect to the argument R, and between derivatives with respect to the argument 
(jl and derivatives with respect to the argument /i/A 2 . 

The new correction to the Hamiltonian gives a contribution of order e 2 to the 
perpendicular and parallel motion of the gyrocenter, comparable to the contribution 

(2) (2) 

from , the term that is usually kept. The correction to the Hamiltonian B gives a 
negligible contribution to the perpendicular drift, but is needed for the parallel motion. 
Thus, both corrections must be kept to obtain the equations of motion to order e 2 . 

It is worth mentioning that equation f)143|) contains the Banos drift |43] in the 
definition of u. To make it clear, instead of choosing the first order Lagrangian as in 
with given in f|lU4j) . we can choose it to have 

T^' = /iV R e 2 ■ ex + ^b ■ V R x b = + /ib ■ V R x b. (147) 



This choice gives a different parallel velocity v! = u — e//b • Vr x b and a different first 
order Hamiltonian if ^ = H ^ + ufib ■ Vr x b. With this new choice, the equation for 
R to first order is 

R=(u' + e/ib ■ V R x b)b + Ax V R B + -^-b x (b • V R b) 

-^V (R±/Ae) (0)xb + O(6 2 ). (148) 

Note that the Banos drift has been made explicit. From here on, we work only with our 
equations of motion f)143j) . (11441) . (I145j) and (1146j) that are equivalent to the equations 
obtained with this alternative choice that makes the Banos drift manifest. 

The gyrokinetic Vlasov equation is readily written for the phase-space distribution 
F(R, u, fi, 9, t) in gyrokinetic coordinates, giving 
dF ■ dF -dF 

_ + R . Vr F + m _ + () _.0, (149) 
or employing the Poisson bracket, 

^ + {F,H} = 0. (150) 

In the absence of collisions and making use of the fact that R, u and 9 are independent 
of gyrophase, the gyrophase independent piece of the distribution function (F) and the 
gyrophase dependent piece F = F— (F) are determined by two decoupled equations |23j, 
namely, 

W + r. VrW+ ^0 (151) 



and 



dt du 



dF dF -dF 

_ + R.V R F + ^ + % = 0. (152) 



Since 9 ~ 1/e ^> d/dt, R ■ Vr, u(d/du), it is usually assumed that any dependence on 
the gyrophase disappears in a very short time and F can be neglected. As a result, we 
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assume that F = (F). By zeroing the gyrophase-dependent piece of the distribution 
function we are eliminating the gyrofrequency time scales, a crucial step to obtain the 
gyrokinetic equation. In the presence of collisions it is possible to estimate the size of 
F and argue that it is small [T7J HI] . 

It is possible to write the gyrokinetic equation in conservative form. Indeed, 
the equations of motion obtained from a phase-space Lagrangian conserve phase- 
space volume. Therefore, the determinant of the Jacobian matrix of the gyrokinetic 



transformation, Jt c , found in Appendix F to be det(JrJ = BT>, satisfies the condition 



7r ' (*n*) + L ^ + le = °" (153) 



For completeness, we prove this equation in | Appendix G| Since 9 and B*^ do not depend 
on 9, this equation reduces to 

Vr • (Bl\R) + §- u (Bfit*) = 0. (154) 

Using this expression and equation (11511) . and taking into account that our choice of 
Lagrangian (11361) implies that dB*/dt = 0, we find the Vlasov equation in conservative 
form, i.e., 

| (B;(F)) + V R • (Btfl(F)) + ^ (B^u(F)) = 0. (155) 
5. Gyrokinetic Poisson's equation 

In Sections [3] and H] we have obtained the gyrokinetic equations of motion and 
the gyrokinetic Vlasov equation keeping the electrostatic potential as an unspecified 
function. The system of equations of electrostatic gyrokinetic theory is closed by 
coupling the gyrokinetic Vlasov equation to Poisson's equation, which is the subject 
of study of this section. 

First, in subsection 15.11 we present the normalization that we employ for Poisson's 
equation. Since several species enter in Poisson's equation, the normalization given in 
Section [2] must be modified. With this new normalization, we obtain the gyrokinetic 
Poisson's equation in subsection 15.21 by simply changing from the coordinates r and v 
to their gyrokinetic counterparts R, u, \i and 9. In Section E] we provide another way 
of obtaining Poisson's equation that is based on field theory. 

5.1. Normalized Lagrangian in a system with several species 
Poisson's equation in Gaussian units reads 

VV(r,t) = -4vr e J d 3 Vi f t (r, v h t) + ^ Z T e J dV/j(r,vj,*) 

-e J d 3 v e f e (r,v e ,t) , (156) 
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where /j(r, Vj,t) is the particle distribution of the dominant ions, Zje and //(r, v/,i) 
are the charge and the distribution function of the impurity /, and f e (r,v e ,t) is the 
distribution function of the electrons. Using the species-independent normalization 

c s t r , A ecz? ,„ „„ s 

t = — ,r=-,A = — — ( 157 ) 

for time, space, vector potential and electrostatic potential, and the species-dependent 
normalization 

v, = ^, V/ = ^, v e = ^, /, = fj = $Lf h U = ^fe (158) 

— v u — v t i — v te — n e0 — n e0 — n e0 



for the velocities and the distribution functions, we get 

- J d^fjv^t). (159) 
Here, n e Q is a characteristic equilibrium value for the electron density and 



Xn e =J T ^ L - (160) 
V 47re 2 n e0 

is the electron Debye length. Recall that e s = p s j L. This new normalization is somewhat 
different from the normalization in Section [2j We have that r = r , A = A, (p = and 
v p = v p (the subindex p can take the value i, e, and any of the values of /), but 
the normalization for the time is different, giving t = rt. The gyrokinetic variables 
obtained in the previous sections can be employed here without further changes because 
both the normalization used in Section [2] and the normalization in equations (11571) 
and (11581) give the same normalization for R p , u p , fi p and 8 p . This normalization is 
species-dependent and for this reason we use the subscripts p = i, I, e to distinguish 
the gyrokinetic variables for ions, impurities and electrons. In the results of Sections [2] 
and [31 the constants e, A, A and r depend on the species, and it will be useful to use 
subscripts for them as well. Their values for ions, the impurity species / and electrons 
are given in Table 15.11 Note the sign in the definition of e e and A e . In Table 15.11 T i0 
and Tio are the characteristic temperatures of the ion species and the impurity species 
J, respectively, and Mj is the mass of the impurity species I. 

The function in ( 170]) is now species-dependent and given by 

4>(R P± / X p e p , R p \\, fi p / \l, 9 P + 7rO(-A p ),t/r p ) 

= 0(R p ±/e s , Rp^fMp/X 2 p , 6 P + 7r9(-Ap), t). (161) 
Later on we will use the the more compact notation 

(p p (R p ,H P ,9 p ,t) = 0( R P ± / e s , Rp\ | , Hp/ Ap, 9 P + irQ (- X p ) , t) . (162) 

To write the total Lagrangian of the system, it is necessary to normalize all the 
particle Lagrangians and Hamiltonians by the same quantity. We choose that the 
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Table 1. Values of the species-dependent parameters e A, A and r for ions, the 
impurity species / and electrons. 

Ions Impurity / Electrons 

Ci = e s /\j e/ = e s /Xi -e e = - e s /A e < e s 

A 2 = \ fT eQ /T ia A/ = Z I ^/T ea m i /T I0 M I -X e = y/mi/m e > 1 

n = ^/T lQ /T e0 ti = v / T I0 m l /T e0 M I r e = ^f mj/m e > 1 

Ai = (T e0 /T 40 ) 3 / 2 A/ = ^Ta/Tmffty/mi/Mi A e = v/^T^T » 1 



normalized Lagrangian and Hamiltonian be £ = C/T e o and H_ = H/T e Q. The gyrokinetic 
Lagrangian (I136p . dependent on the species, becomes in the new normalization 

dRp T p0 e s d6 

-Vp-jr ~ h p ( 163 ) 



T, 



P rp 

J- eOTp 

where 



^A(R P ) + u p b(R p ) - ^ P K(R P ) 



dt T e ()XpTp — dt 



T, 



FT - ~ p0 
H P ~ 



^ Mp 2 + fi p B(R p ) 



T e 2 

+ ( r p> «p> fh) + Z P e s(4>)(R v ±/e s ,Rp\\,iJ l p/\l,t) 

1 P nA 



-eO' 



+ 



Z n A„e 2 



p^pt s lTf (2) 



A, 



^ > {R p± /e s ,Rp,^p/X 2 p ,t_, X p ) + (R p± / 6 „ R p , u p , /y /i p /Ag, t, X p ) (164) 



and 



K(Rp) = -b(Rp)b(Rp) • Vr^ x b(iy - V^e^Rp) • e.iRp). (165) 



Here, obviously, Z{ — 1 and Z e = —1. Analogously to (I162p we define 



^^(R-p'/V^) — ^ ( <h' (R P ±/ts,R P , ^p/ X p ,t, X p ) 



,(2) 



(166) 



and 



^l^i^^^t) = V%(Rp±/e s , Rp, Up, /ip, ij, p / Ap, t, X p ) . (167) 

The equations of motion are those obtained in (I143p . (11441) . (I145p and (11461) 
multiplied by r p . As a result, the Vlasov equation for each species is 



dF p . dF p 
^ + t p R.Vk^ + V^==0. 



:i68i 



In what follows we work in non-dimensional variables but do not underline them. 



5.2. Gyrokinetic Poisson's equation via the gyrokinetic change of coordinates 

Our objective is to write Poisson's equation (I159p in terms of F P (R P , u p , /i p , t), i.e. the 
solution of equation (I168p . Sinc^3 (r, v p ) = T ep P (R p , u p , fi p , 9 P , t), 

Fp(Rp, Up, fip,t) fp(T tpt p(Rp, Up, fip, 9p,t) ,t) . (169) 

* We stress that the transformation T e depends on the species through the values of e p , A p and A p and 
denote it by T CpiP . 
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(170) 



= j d 3 Riduidnid6i |det (Jr e<ii ) \ Fi ^(V (r eij j(Rj, u h fi h U t) 

+ ' S ^Z[ j d 3 Riduidfiid9j det ^ Jt e/ 7 j Fi 8\7r r \T ei j(Hi,ui, fii,9i,t) 
i 

- j SRAuA^M, |det (J T ,J | F c 6 (r„(H., «., „«, f) 



(171) 

where Jt^ p>p (R p ,u p , fi p ) is the Jacobian matrix of the transformation T £piP , calculated in 



Appendix F and 7r r (r, v p ) := r is the projection onto the spatial part of the coordinates 



(r, v p ). From an abstract viewpoint this is, perhaps, the simplest way of writing (I159p in 
terms of F p (R p , u p , /i p , t). Observe that (I17ip is an exact relation. However, in practical 
terms, we have only computed the explicit expression of T ep)P up to a certain order in 
e p , or equivalently, in the species-independent parameter e s . Making use of the results 
of this paper we can give an explicit expression for the change of variables up to order 
e 2 , namely 



• I F epj p(Rp, Up, /ip, 8p, t) 

,2 



Rp + — p(Rp, /ip, 9 P ) 

An 



e 

+— 



r> , „ d P ,a d P 

Odp A 



A dfip 

which allows us to write the identity 
S ( T ^,p( R p' U P' /V ^P' *) 



+ 0(e 



3 ) 

s I • 



(172) 



+ 



= 5 (tip + -^p(Rp, /ip, P ) — 

^(R P ,2 + ^,i^- + ^,i^ ■V5(Rp+^p(Rp,fip,0p)-^+O(e 3 s ), (173) 

where V<5 denotes the gradient of the Dirac delta function with respect to its natural 
arguments and the subscript p = i, I, e in the corrections Rp,2, A*p,2 and Pt 2 indicates 
that these corrections depend on the species through A p , X p and P (R P , fi p , P , t) in (I162p . 
Substituting (I173p into (11711) finally gives 



L 2 



VV(r,i) = Z~2 Z p d 3 RpdupdfipdO p B^ p Fp 6(r p + ^p(Rp, /x p , P 

p LP 

dp 



+ . . . (174) 



Here we have used Appendix F| to write det( Jt^ p , p ) = B\\ p - We employ the subscripts 
p = i, I, e in the determinant of the Jacobian B! because it depends on the species 
through its dependence on e p = e s /\ p . Expression (11741) may seem accurate to order e 2 , 
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V8[R p + fp-r) = [ I + ^V Rp p) .V Rp 5(R p + ^p-r 



i.e., that on the right side of ( I174|) we are dropping only terms which are of order or 
higher. However, it is easy to see that this is not true. Using 



(175) 
(176) 



A 



A, 



V Rp 5 \R p + ^-p-rj +0(e s ), 



and 



e s dp 

e s dp 

X p dOp 



VS ( R p + -^p - r 

A p 



V<5 ( Rp + y-p - r 

An 



_d_ 

dfip 



d_ 



5 ( R p + — p - r 

Ap 



5 ( Rp + — p - r 



(177) 



and integrating by parts so that the delta function does not appear differentiated, we 
find that (I174p becomes 



e sX 2 De 
L 2 



V V(r, t) = Y, z P J d 3 R p dupdfipd9 p 6 (r p + ^p(R p , fi p , 6 P ) - r) j 

2 ) + ^K^0 



^vp |_/»p 

+ wS B * ] ' pFp6pA 



178) 



Note that terms that seemed to be of second order in e s are in reality first order 
contributions. Similarly, terms that seem to be of third order and are neglected in 
f ll74p are in reality of second order. To obtain a gyrokinetic Poisson's equation correct 
to order e 2 it is necessary to carry (11721) to an order higher 

ii i 7 1 £;)j p(Rp, tip, \Xpi Qp-, l) 



Rp + — p 



+ ( Rf>,3 + Rp,2 ' V Rp P + fJ, p> 2 



d 2 p 1 2 d 2 p 



^R P ,2 - 


dp 




dp 

2 n + 
dfip 




i 2 

2^^ 


+ 0(£ 







(179) 



d^ipdOp 

All the terms entering this equation are computable from the results found in Section [3J 
We leave this for future work. 

By integrating over the delta function in (I178p we make the gyrokinetic Poisson's 
equation more explicit, 



e sX 2 De 
L 2 



VV(r,t) Z p f dupdfXpdOp ^ - 

p J det(I + (e s/ 



/Ap)V Rp p) 



B Ip F p 



d 



d 



wA B \\« Fp9p > 1 



dfj, p 



B\\, p Fpf J 'P,i 



(180) 



Rp— R(r ,p p ,9 P ,e s / X p ) 
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where R(r, fi p , 9 P , e s / X p ) is the the function defined by solving for R p the equation 
R P + {e s /X p )p(R p ,fJ, p ,9 p ) = r. 

Finally, we would like to point out that usually equation (1180p can be simplified 
even more. The following discussion is not meant to be an exhaustive review of all 
possible orderings and simplifications, but a brief comment on the most typical approach. 
Usually F p has contributions with different characteristic scales of variation, ranging 
from the shortest scale p s to the background profile variation scale L. At the same 
time Vr p Fp ~ 1 because the saturation amplitude of the short scale fluctuation F PtPs 
is sufficiently small, i.e., F P:Pa = e s F PtPs (R p± /e s , R p \\,u p , fi p ,t) ~ e s <p(r±/e s ,t), where 
F P)Ps is a function of order unity with derivatives with respect to its arguments of order 
unity. For Vr p F p ~ 1, we can use the approximate expression R(r, /z p , 9 p , e s /X p ) = 
r - (e s /A p )p(r,/i p ,6> p ) + 0(ef) because F p (R(r, \i p , 9 P , e s /X p ), u p , \i p , 9 P , t) = F p (r - 
(e s /X p )p(r, n P ,9 p ),u p , (j, p ,9 p ,t) + 0(e 2 ). In addition, the term containing R Pi2 in (11801) 
becomes of next order because Vr p F p ~ Vr p -B|J p ~ Vr p • R Pi 2 ~ 1. To simplify even 

more, we use det(I +(e s /A p )V Rp p) = 1 + (e s /A p )V Rp ■ p + 0(e 2 ). As a result, 
CsA ^VV(r,t) = J2 Z pJ du p dfx p d9 p \Bl p F p - ^ Bl p F p (v Rp • p) 



V- 

I N 



(B* ]p F P n Ptl ^j + — (B* ]p F p 9 Ptl ^ |^ 



+ 0(e 2 s ). (181) 

-(e a /A p )p(r,/ip,e p ) 



dfj, P 

This is the simplest Poisson's equation correct to first order in e s . As we have already 
pointed out, it is necessary to include the corrections R P) 3, [i p ^ and 9 P)2 to find the 
contributions of order e 2 s . We finish by checking that the assumption Vr p -F p ~ 1 is 
consistent with the gyrokinetic system of equations. Consider the short scale pieces of 
the Vlasov equation (11681) and of the gyrokinetic Poisson's equation (I18ip . The result 
will be a typical Sf formulation [21 IH El El E] , demonstrating that it is possible to find 
a closed non-linear system of equations to determine the short scale fluctuations F PiPs 
and if. This does not mean that Sf formulations are always valid, but it indicates that 
the assumption Vr F p ~ 1 is consistent with the gyrokinetic system of equations. 



6. Gyrokinetic Field Theory 

Gyrokinetic field theory is the formulation of gyrokinetics as a classical field theory and is 
defined by an action functional S that depends on the trajectories of the particles and 
the electromagnetic field. A closed system of equations coupling the electromagnetic 
field and the gyrokinetic distribution function is obtained by finding the stationary 
points of S. The first application of field theory to plasma physics is the work by Low 
in jl5] that was later extended to gyrokinetic theory in the seminal papers by Sugama 
and Brizard [17] ■ The development of gyrokinetic field theory in the last decade 
has been motivated mostly by the fact that it allows to identify in a systematic way 
conservation laws from symmetries of the Lagrangian. For the sake of completeness 
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we briefly recast now our results in field theory language following the Lagrangian 
formulation of Sugama [46J. 

In subsection 16.11 we show how to obtain the gyrokinetic equations of motion and 
the gyrokinetic Poisson's equation from a variational principle. The advantage of this 
procedure is that it ensures that an energy-like invariant is conserved. We obtain the 
invariant and prove that it is conserved in subsection 16.21 



6.1. Equations via variational principle 

Using the abbreviated notation Z p = {R p , u p , fi p , 9 P }, the action functional is 



: [z p (., 

+ J2[ 1(lt / ^^o%(Z p o)i ? P Q(Z p o)^(z p (Z p0 ,to;*),Z p (Z p o,tQ;t),^(«,t)), (182) 



where F p0 is the distribution function of species p at time t , and Z p (Z p0 , t ; t) is 
the trajectory in phase-space of a particle of species p satisfying the initial condition 
Z p (Z p0 ,t ;t ) = Z p0 at t = t . The first term on the right side of f l 1 8 2 j) is the action 
for the electric field (recall that we are considering a static magnetic field) and the 
second term is the sum of the actions of the particles. The gyrokinetic Lagrangian C p 
of the species p is defined in (I163p (recall that we have dropped the underlining for 
normalized variables). Its third argument stresses that C p depends as a functional on 
the electrostatic potential. 

The Euler-Lagrange equations for £ are obtained by finding its stationary points 
under infinitesimal variations of the maps Z p (Z p0 , t&t) and ip(r,t). The allowed 
infinitesimal perturbations to y?(r,t) vanish at the boundary of the spatial domain of 
interest, and the perturbations to both Z p (Z p0 ,to;t) and (p(r,t) must be zero at t — to 
and t = t\. The calculation of the variation with respect to Z p (Z p0 , t ;t) gives the 
gyrokinetic equations of motion and is a repetition of that leading to equations (I143p . 
(11441) . (I145p and (I146p (recall the discussion on the phase-space Lagrangian methodology 
in subsection 13. ip . The distribution function at time t is 

-fp(Z P , i) := F p0 (Z p o(Z p ,t;t )), (183) 

where Z p0 (Z p , t; t ) is the inverse of the map Z p (Z p0 , t ] t), i.e., Z p0 (Z p (Z p0 , t ; t), t; t ) = 
Z p q. Note that F p automatically satisfies (11681) . Observe also that from condition (I153p . 
which is automatically satisfied by (I143p . (I144p . (I145P and (I146p . we obtain 

-B|* )P (Z p o(Zp, t; t )) 

|det(J Zp0 ^z p (Zpo(Z p ,t;t )))| 

where (Jz p0 ^z p )a = ^Z^/dZ p0 is the Jacobian matrix of the map Z p (Z p0 , t ;t), and 
det(J Zp0 ^ Zp (Z p0 )) = £|* jP (Z p0 ) / 'B^ p (Z p (Z p0 ,t ;t)) is its determinant. This property is 



d 6 Z p0 Bl v (Z p0 ) = d 6 Z p TJ - zrF ^ P frj P \„ -Tl = ^^^||,p(Zp), (184) 



proven in Appendix H for completeness. 
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The stationary points of S under variations of <p(r,t) are given by S^Yl = 0, with 
5 V E = J dtd 3 rV5<p(r,t) ■ V<p(r,t) 

— ^^Z p e s / dt d 6 Z p0 B^ p (Z p0 )F pQ (Z p0 )5 lfi ((f) p (Z p (Z p0 , £ ; t), t)) 

v J 

— ^ P ^ P s / dtd 6 Z pQ B* p (Z p0 )F p0 (Z p o)5 v ^^ p (^Z p (Z p0 ,to] t),tj 

p p 

7 2 p 

~z2^r / ^^^o5i|, p (Z p o)^o(Z p o)^^l P (z p (Z p0 ,to;t),t). (185) 

v v 

In Appendix I we have evaluated all the terms. The final result is 

8 V H = e s J dt d 3 r 5tp(v, t)V{r, t), (186) 

where 

?>M) = -%?vVM) 

~Z~2i Z p J d 3 Rpdu p dfi p d9 p B^ p (R p , u p , /i p )F p (R p , u p , fi p , t)^S (r p + j-p - r^j 



^-R p ,2± • V Rp 5 ^R p + j-p - r^j + /-W^-<* (ftp + Y p 



A p 

+ e ^W p 6 fa + T/-* 



187) 



By imposing that 8 V H = for any S(p(r,t), we find that V(r,t) = 0. This is Poisson's 
equation. By integrating by parts to leave the delta function undifferentiated, and 
then integrating over the delta function, we recover (11801) . Note that even though the 
Hamiltonian is obtained to second order in e s , and hence the Vlasov equation is also 
known to second order, the gyrokinetic Poisson's equation that we have found using the 
variational principle is only correct to first order in e s . It coincides to first order with 
the first order equation that we found independently in subsection 15.21 

Observe that equation (I187p keeps some second order terms that could have been 
neglected, as done in (11811) . These terms are important if we want to take advantage 
of the field theory formulation of gyrokinetics. The action (I182p is invariant under 
time translations, so Noether's theorem automatically provides a conserved quantity, 
interpreted as the total energy of the system, which we denote by rl(t) because it is 
indeed the field theory Hamiltonian. However, 'Hit) =0 on the equations of motion of 
E, that is, on solutions of (jT33]» . fTM . ([USD , flBBfl and ffTSTj) . without neglecting any 
terms. We show this in the next subsection. 



6.2. Conservation of energy 

In this subsection we prove that if the equations of motion are obtained via a variational 
principle of the action ()182p . there is an energy-like invariant T-L(t). This is an application 
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of Noether's theorem, and it is based on the fact that in the action (I182p the only time 
dependence is through the functions Z p (Z p0 , t ; t) and (p(r,t). To obtain the energy- 
like invariant, we use the equations of motion derived from the variational principle in 
subsection 16.11 

In order to find the conserved quantity we first perform the change of variable 
t — t' + St in every integral on the right side of (11821) . which, of course, does not change 
the value of the action, giving 



tx-St 



dt' 



t -St 



\ 2 p 2 
A De e s 

2L 2 



d 3 r \V<p(r,t' + 8t)\ : 



+ J d 6 Z p oB*^ p (Z p0 )F p0 ( y Z pQ ) C p \ t=t ,_ 



■st 



where 

z 



P\t=t'+St 



Cpi Z p (Zp , t ; t' + St), Z p (Z p0 , t Q ; t' + 5t), tp(; t' + St)\. 



189) 



We expand the right side of the previous equation up to first order in St to find 



"TV > ! /) 



+ <L E 



-TV > i )i 



<P(* 



S(p=5t(d<p/dt) 



-st 



r x 2 ^ 2 

A De e s 
2L 2 



d 3 r\Vp(r,t')\ 2 + J2 / rf 6 ^ P o5||, p (Z p o)^o(Z p o)^| t=t , 



n *'=ti 



(190) 



t'=to 



where <^,E| stp=st(d<p/dt) an d ^ZpE^z are the variations of E under perturbations of 
both <£>(r,t) an d Zp(Z p o,^o;^) with the specific form Stp = St(dip/dt) and SZ P = StZ p . 
For the variations with respect to Z p , we obtain 

^z.,E = / dt' d 6 Z p0 B* p (Zi p q)F p0 (Zi p q)5z p C p (^Zi p (Zi p q, t ; t'), Z p (Z p0 , t ; t')j 



to 



to 



<tf / rf 6 Z p0 5 | ] jP (Z p0 )F p0 (Z p0 )^ ( + 

a=l 



dZ 

d 6 Z p0 5|* p (Z p0 )F p0 (Z p0 ) 2^ 7T^~ 

«=i " Z p 



t=t' 



v 

i t'=ti 
t'=t 



dZ« 



t=t' 



(191) 



where to obtain this last equality we have integrated by parts in t and we have used the 
equations of motion 



d I dC n 



dt \ dZ°t 



dZ° 



a = 1, 2, . . . , 6. 



(192) 



Applying SZ p = StZ p we find 



6*JE\ 



St 



■ dZ 

d G Z p0 5|| p (Z p0 )F p o(Z p0 ) 22 Z p -^rr 



1 tf=ti 



t=t' 



(193) 



t'=t 
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For the variations with respect to Sip we obtain that 

where we have employed (I186p and Poisson's equation V(r, t) = 0. To obtain (I186p we 
assumed that Sip vanishes at the boundaries of the domain, and this is not necessarily 
the case for Sip = St(dip/dt). We avoid this problem by assuming either that the domain 
extends to infinity, where ip — 0, or that we are in a periodic box and the contribution 
from one half of the boundary cancels with the contribution from the other half. In 
either case, we are assuming that there is no net energy flux through the boundary. 

By substituting (I193p and (11941) into (11901) and observing that the term linear in 
St in (I190p has to be identically zero, that is, the sum of the second, third and fourth 
terms in the right side must vanish, we find 



St 



A De e s 
2L 2 



d 3 r \V(p(r,tf) 



d e Z p0 B* p (Zp )Fpo(Z p0 ) 



h dz ? 



t=f 



0. 



Then, 



U(t>) 



1 *'=*! 



0. 



(195) 



(196) 



for solutions of the equations of motion, where 



u(t) = j2 J d * z p» B \\ 



po d \\, p K^po)F v q{'Za p q)H p 



Z p (Z p0 ,t ;t),p(«,t) 



A 



De c s 



2L 2 



J d 3 r\Vip(r,t)\ 2 . 



Since this must hold for every interval [to^i], we deduce that 
dU{t) 



dt 







(197) 



(198) 



on the equations of motion. Finally, using equations (I183p and (I184p . we can write the 
integrals over Z p0 in (I197p as integrals over Z p , giving 

H(t) =J2 I d 6 Z p Bl p (Z p )F p (Z p ,t)H p (z p ,ip(.,t)) - ^ I d 3 r\Vip(r,t)\ 2 . (199) 

p J J 
The conservation of this energy-like invariant is only satisfied on the equations of motion, 
that is, the variational equations of motion and the variational Poisson's equation 
obtained in subsection 16.11 must be used. It is important to keep all the given terms, 
even if they are higher order than first (the order to which Poisson's equation is correct). 
For example, it is necessary to obtain the exact function R(r, /i p , 9 P , e s /\ p ) in (11801) to 
have exact conservation of "H. 
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7. Conclusions and further work 

The gyrokinetic ordering in a static magnetic field consists of the ordering assumptions 
(jTJ), defined by a single parameter e = p/L = uj/Vt. In this paper we have strictly 
implemented the gyrokinetic ordering in the phase-space Lagrangian to obtain explicitly 
the gyrokinetic Lagrangian to order e 2 for general magnetic geometry. 

Our approach differs from previous phase-space Lagrangian (or Hamiltonian) 
derivations of gyrokinetics. In previous work [25] the calculation is performed in two 
steps. First, with zero fluctuating electrostatic potential, an expansion in powers 
of e = p/L is performed and a gyrophase independent guiding-center Lagrangian 
is determined to order e. Then, the electrostatic fluctuations, whose size is given 
by a new expansion parameter e v = Zetp/Mvf, are switched-on, reintroducing a 
gyrophase dependence that is removed order by order in yielding the final gyrokinetic 
Lagrangian, usually computed up to order e 2 . When the expansion in e v is performed, 
the fact that there has been a previous expansion in e is ignored and the terms of order 
ee v are never calculated. Thus, the final Lagrangian is missing relevant terms of order 
ee v and e 2 . 

The novelty of our work can be easily understood by examining the explicit 
expression of the second-order gyrokinetic Hamiltonian, H , given in equations (I132p . 
fll33p . fll34p and (I135p . It shows in a transparent way that gyrokinetic theory ties 
together geometry and turbulence, so that no splitting between guiding-center and 
gyrokinetic dynamics is possible. See, for example, in (I134p where magnetic 

geometry and electrostatic potential appear together. This is the first time that the 
electrostatic gyrokinetic equations in general geometry are fully computed to order e 2 

(2) 

and the calculations are pursued to the point of reaching formulae, like the one for H , 
that can be straightforwardly implemented in a computer code. 

From the new phase-space Lagrangian in (I136p and the new Hamiltonian in (11371) 
we obtain a new Vlasov equation and a new gyrokinetic Poisson's equation. The Vlasov 
equation is correct to second order in the expansion parameter e, and it is to our 
knowledge the highest order full / gyrokinetic equation available in the literature for 
general geometry. In the limit where the electrostatic potential has a scale of variation 
much larger than the gyroradius of the species of interest, this equation is also the 
highest order drift kinetic equation that we are aware of. 

The gyrokinetic Poisson's equation derived from the new Lagrangian is, however, 
only correct to first order in e. We have calculated Poisson's equation employing two 
methods. In Section [5l in the integrals in velocity space that enter in Poisson's equation 
we have simply changed from the coordinates {r, v} to the gyrokinetic variables. In 
Section [H] we have used a variational formalism. The equation obtained with the 
variational formalism is set by the form of the Hamiltonian and it is only correct to 
first order, although it contains terms that are higher order. These higher order terms 
do not add accuracy, but they are necessary to have an exact energy-like invariant. By 
directly changing from {r, v} to the gyrokinetic variables in the integrals that enter 
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Poisson's equation, it is possible to obtain a higher order Poisson's equation if the 
higher order corrections R 3 , /i 2 an d 2 are calculated. It is not necessary to calculate 
the Hamiltonian to next order. The calculation of R3, fa and 62 is however very tedious 
and is left for future work. If this procedure is followed and the third order Hamiltonian 
is not obtained, the conservation of the energy-like invariant will not be exact. 

Some natural steps following the present work are the extension to the 
electromagnetic case and the introduction of external flows. We will also investigate the 
implications of the new terms in the Hamiltonian for the transport of toroidal angular 
momentum in tokamaks. 
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Appendix A. Proof that the inverse of the matrix in (120ft defines a Poisson 
bracket 

Let L a/ g(Z) be an invertible, antisymmetric matrix of dimension 2n whose components 
are functions defined on a region U C M 2n , and such that 

Denote by P^(Z) the inverse of L a0 (Z), P a ?(Z) = (L^Z))^. Then, P a P defines a 
Poisson bracket by contraction with the differentials of pairs of functions on U, i.e., 

We must prove that (jA.20 satisfies skew-symmetry, given in (1241) . the Leibniz rule, 
given in (T25]) . and the Jacobi identity, given in fT26|) . Skew-symmetry is satisfied because 
the inverse of an antisymmetric matrix is antisymmetric, i.e., P a/3 = —P^ a . The Leibniz 
rule is trivial to check. As for the Jacobi identity, it is immediate to see that it is 
equivalent to 

£^w + £^w + I> 7 w = < '^''^ k (A ' 3) 

5=1 8=1 5=1 
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Showing that flA.3j) is equivalent to flA.lj) is a simple exercise of application of the 



formula of the derivative of the inverse of a matrix, 

^!! = _f;( L -. r ^( L - r . (a.4) 

5,p=l 

Appendix B. Calculation of the Lagrangian after the non-perturbative 
change of variables 



In this Appendix we prove equations f l83|) and f l84j) . For equation (T83T) . we use that 
according to flHTT) 

Vr^SW = - / /f [VA(R, + ep') • p' + eV Rs p' ■ VA(R, + ep') ■ p' 

+V Rg p'-A(R g + ep')}. (B.l) 
Recall that the prime ' here indicates that the function depends on fi' g . Employing 

1 T VA(R g + ep') ■ p' = -±jp' • VA(R g + ep') + -Lp' x [V x A(R, + ep')} 



2n' g " y r ' 2p' g r ^ r ' 2^ 

X c) 1 

- T [A(R g + ep')] + —p'xB(R g + ep'), (B.2) 



2^ 



e^ 1 v 9 ^ n 2^ 
e 7 V R9 p'-VA(R 9 + ep')-p / 



^7V Rs p' • [p' ■ VA(R g + ep')} + ^V Rff p' • [p' X B(R 3 + ep')} 
V Rg p' • ^[A(R, + ep')} + ^7Vr,p' ■ [p' x B(R, + ep')} (B.3) 



and 



L VrsP ' ■ A(R 9 + ep') = V Rff • A(R g + ep'), (B.4) 



equation (IB.lj) becomes 



Vr 9 5jvp = — -A(R g + ep) + ^A g - V Rg p ■ A(R g + ep) 

[p' x B(R, + ep') + eV Rff p' ■ (p' x B(R 3 + ep'))] . (B.5) 

'0 Z ^g 

To obtain equation (153"]) from this equation we use (1751) to write 
da' 

ffp'xB s = pxB 9 , (B.6) 

'0 Z H'g 

and we employ (177)) to get 

V Rfl p' ■ (p' x B(R g + ep')) = [(p' x b 9 ) • B(R 9 + ep')] V Rs b g • p' 

+-^[b 9 • B(R, + ep')]V Rg e 2g • e lg . (B.7) 
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For equation (1841) . we use that according to (JHT 



OS 



NP 



I'g 



d9 9 
Employing 
e dp' 



2// 89 9 



W 9 

VA(R 9 + ep') ■ p' 



e^-VA(R s + 6 P ')-p / +|^-A(R g 



f ° pl [p'-VA(R g + ep')} + ^-^-[p'xB(R g + ep')] 



and 



^^•A(R 



•A(R ff + ep'), 



(B.8) 



(B.9) 



(B.10) 



• a ~« a d9 g dfi' g " K ~ 9 

equation flB.8|) becomes equation (jgljh To obtain the final form of the equation we have 
also used (1751) to write 



^•[p'xB(R 9 + f p')] 



VB(R s + ep' 



(B.11) 



Appendix C. Comparison to first order with the results in reference | |17| 

In this Appendix we compare the gyrokinetic variable transformation obtained in this 
article with the variables found in [T7] . 

To be able to compare with the results in [T7J, given in the form Rpc(r, v, t), 
Epc(r,v,t), /xpc(r, v,t) and 9pc{r,v,t), we will use the transformation (r, v) = 
T e (H,u, fj,,9,t) to write them as Rp<?(R, u, /i, t), E PC (R,u, fi,t), fipc(R,u, fi,t) and 
9pc(R,u, fi,9,t). After doing so, we will see that the gyrokinetic variables Rpc, Epc 
and fipc in [T7| are gyrophase independent quantities up to the order that they are 
defined. This is a property that must be satisfied because otherwise the new gyrokinetic 
variables would have fast time dependence through the gyrophase. The variable [ipc 
must be a function of /i only because there is only one adiabatic invariant associated 
with the gyromotion. 

The gyrokinetic variables of [TTJ are the gyrocenter position 

2 



R P (7 = r H v x b H 

B B 



ti||b + j (v x b) + (v x b) (v\\h + 



B 



e 2 v\\ 



: -^-Vi ■ Vb + -^bb • Vb ■ v ± + ^-^b[v ± v ± - (v x b)(v x b)] : Vb 

Ae 2 



V (R± /Ae)$ x b + 0(e 3 ), 



A5 2 

where abi M=ax (b- M); the gyrokinetic kinetic energy 
£ PC = ^ + Ae0 + O(e 2 ); 



(C.l) 



(C.2) 
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the magnetic moment 

^=i"||(vx6)-V B -||b.Vb.(vxb)-^b.Vxb 

' ''" [yi(v x b) + (v x b)vj : Vb> + ^0 + 0(e 2 ); (C.3) 



AB 2 ' K ' ^ ' B 

and the gyrophase 

/v • e 2 \ e euji - e 

6» PC = arctan — - t^v_l ■ V5 ^-b ■ Vb ■ v± + — (v x b) ■ Ve 2 ■ ei 

\v • eij B 2 v\B B 

" " v ± v ± - (v x b)(v x b>)] : Vb - Ae ®® + 0(e 2 ). (C.4) 



Av\B v x ± v A ;J A 2 5<9(/i/A 2 ) 
To these variables we apply the gyrokinetic transformation. First we apply the 
non-perturbative change of variables T/vp e . For the variable Rpc we find 

R PC = R 9 -eV-V Rs Mx (v l{g b g + p x B 9 ) + x b a ) ■ V Rfi b a 



6 2 



= V Ro 



,-2„,.. 9 



+ ^bb • V Rs b ■ (p x b p ) - ^b g [pp - (p x b 9 )(p x h g )\ : v Rff b, 



-"9 

Ae 2 „ . , .,, 



\Bl 



V (R±/Ae) $xb g + 0(e 3 ), (C.5) 



<-) 



where we have used r = R 9 + ep, B(r) = B g + 0(e), b(r) = h g + 0(e), i>|| = + 0(e) 
v_l = p x B g + 0(e) and 

v x b(r) = i-v x b, - ep ■ V Rg I ^ I x (i^b, + p x B ff ) + 0(e 2 ) 



-p-ep-V Rs (M x ( W || 9 b 5 + pxB 9 )+0(e 2 ). (C.6) 



b 

h 

We then apply the transformation (R fl , v\\ g , fi g , 9 g ) = Tp >e (R, u, //, 9, t) by using H g = 
R+e 2 R 2 +0(e 3 ), with R 2 given in ffTTUl) . = w+0(e), // s = /i+0(e) and 6 g = 6+0(e). 
As a result we find the expression 

2 2 

Rpc = R + ^-b[pp + (p x b)(p x b)] : V R b + ^[pp + (p x b)(p x b)] ■ V R S 

A A A / b \ ^ 

+ — (p x b) • V R b + e 2 wb • V R I - I x p - — bb • V R b ■ (p x b) 
--bxV R b.p + 0( e 3 ), (C.7) 

£5 

where i? and b depend on R. To simplify expression (10. T|) we use 

pp+(pxb)(pxb) = -|(I -bb) (C.8) 
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and 

(p x b) • V R b - b x V R b • p = [(b x V R ) x b] x p = (V R ■ b)(p x b). (C.9) 
We finally obtain 

9 9 

R PC = R - ^ybb • V R B + ^V R± £? + 0(e 3 ). (CIO) 

To find this result we have used Vr • b = — B^b ■ VrB. 

Applying the non-perturbative change of variables to Epc we find 

Epc = \v\ g + \\p x B g | 2 + Ae0 + 0(e 2 ) = l -v\ g + fi g B g + Ae0 + 0(e 2 ). (C.ll) 

To transform to the variables R, u, /i and 6 we use v\\ g = u + eu\ + 0(e 2 ), /x s = 
/i + e/ii + 0(e 2 ), with «i and \L\ given in (II 111) and fl 1 1 2 1) . R 9 = R + 0(e 2 ) and 
g = + O(e), giving 

£pc = ~it 2 + /iS(R) + 0(e 2 ). (C.12) 
For /ipc, we find that the transformation (r, v) = T N p e (R, g ,v\\ g , (j, g ,6 g ) gives 

Vpc = fi g - ev llg p ■ V Rs b g • (p x b 9 ) + -^b 9 • V Rg b 9 • p - f!^b ff • Vr 9 x b g 

+ ^[(p x b 9 )p + p(p x b,)] : V Rs b g + ^0 + 0(e 2 ), (C.13) 
where we have used = v\\ g + 0(e), v_|_ = p x B 3 + 0(e), 

^ r ^. (v . Ul rf\-^ P -V R ,B, + 0^ (CM) 
and 

t; 2 - (v ■ b(r)) 2 = vf lg + 2/i g B g - [(v ]lg b g + p x B ff ) • b(r)] * + 0(e 2 ) 

= 2fi g B g - 2ev llgP ■ V Rg h g • (p x B g ) + 0(e 2 ). (C.15) 

Then employing fi g = \i + e/ii + 0(e 2 ), with /ii given in (11 12j) . R 9 = R + 0(e 2 ), 
v\\ g — u + 0(e) and 9 g = 9 + 0(e), we find 

A*«7 = A* " ^^b, • Vr 9 xb 9 + ^[p(p x b s ) - (p x b 9 )p] : V R9 b 9 + 0(e 2 ) 

= /i + 0(e 2 ), (C.16) 

where to obtain the final equality we have used (1(180 . 

Finally, for 8pc we obtain that the non-perturbative transformation gives 

= 9 g - -^p ■ V R9 b, ■ p - — (p x b g ) ■ V Rg B g - ^-b g • V R9 b 9 ■ (p x b ff ) 

' '>[pp - (p x b g )(p x b,)] : V R9 b 9 - + 0(e 2 ), (C.17) 
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where we have used 



v • e 2 (r 

arctan | — f | = arctan 

v ■ ei(r 



€ 



(v llg b g + px Bg) -e 2 (r) 
(v llg b g + pxB g )- ei(r) 



2 MgJ B, 



■[(p x B g ) ■ e 2g ]p ■ V Rs e l9 ■ {v\\ g h g + pxB 9 ) + 0(e z 



9 
61^1 1 

= d g - -^p ■ V Rs b 9 ■ p + ep • V R9 e 29 • e l9 + 0(e 2 ). (C.18) 
l\i g 

Here we have employed V Rg ei 9 = -(V Rg bg ■ e lg )b g - (V Rg e 29 • eig)e 2 g and V Rs e 2;; = 
— (Vr l3 g • e 2g )b g + (VR g e 2 g • e lg )e lg to write the last equality. Equation (1C17[) is now 
rewritten by using 9 g = 9 + e6»i + 0(e 2 ), with 9 X given in ffTTH]) . R g = R + 0(e 2 ), 
= u + 0(e) and n g = y. + 0(e), giving 

9 PC = 0- —\pp + (p x b)(p x b)] : V R b + 0(e 2 ) = 9 + —b • V R 7? + 0(e 2 ), (C.19) 

where we have employed (1C.8P and Vr • b = — B^b ■ Vr5. 

From equation flCHU]) . (jHl2j) . flCTTBl) and flClT9|) we see that the gyrokinetic 
variables defined in [T7| are simple functions of the gyrokinetic variables in this article. 
Notice that Rpc and Epc do not depend on the gyrophase and that fipc is equal to \x 
to the order of interest, as expected. 

Appendix D. Calculation of the second order gyrokinetic Hamiltonian 

In this Appendix we show how to obtain f 1 1 3 2 j) from equation (11301) . 

First, by employing equations (19T]) . (193]) . fTTU8"j) . fTTTU]) and (17231 . we find that 
(if) = 0, (Vr4 2) ) = 0, (dsf/dt) = 0, 



and 



( F i 2) ) = -bb) : VrVrB • b. (D.2) 



Here we have used equation fl 102[) repeatedly. These results give that equation (I130p is 
equal to 

H i2) = -f^b ■ V R b ■ V R 5 + i^(I -bb) : VrVrB ■ b - ^|Vr^| 2 + <|> 

+ (r 2 ■ v (RxA )F (1) ) + Ti + r 2 + r 3 + r 4 , (d.3) 

where we have used (Vr x b) x b = b • Vr£> to write 

b • [(v R x b) x (r 2 )] = ^[(v R x b) x b] • v R 5 = ^b • v R b ■ v r b. (da) 

The terms 71, %, Tz and 74 in equation (1D.3|) are 



Ti = fi( (r^ + R 2 -V R A) (D.5) 
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T 2 = B 
T 3 = B 



and 

% = - u b ■ ( m R 



R 



b +tti) b 



9R 2 
~~09 



r (1) // 

1 e ~ ^1 



OT 



(i) 

R 



+ 9 



~09 
OT 



45 
(D.6) 

(D.7) 



(i) 

R 



Ou On 09 



+ B(fn 



OT 



(i) 



dfj, 



9r (i) 



(D.8) 



We proceed to write these terms in more convenient forms 
The term Ti in (1D.5j) is rewritten using 

9R 2 1 



R 2 • V R A 



to obtain 



89 ' 2 
1 



R 2 -V R A.^-^-V R A.R 2 



09 



09 



ld_ 
209 



9R 1 
2V BxR 2 ).^ + --(R 2 .V R A.R 2 ) 



r 1 = B 



-a) 

R_L 



1 



B x R s 



9R 2 



Employing (|47p . this last result becomes 



1 



r\ = §( (r«-v ( K, /e) 4 2) 



b X 



— — + V 

90 



95 



(2)' 



CRi/e) 



09 



[TL 2 ■ V R A • R 2 ) 

(D.9) 

(D.10) 
(D.ll) 



Realizing that 



R 



b X V(R ±/£ ) 



9S 



(2)' 



9_ 
90 



09 



V(R_ L / e )S'^ 2) • 



(2) 



u 9T^ 



we finally obtain 



71 



^(i) 

R 



V(R ± / e )S'^ ) 



b x V(R x/e) 



95 



(2)~ 



90 



(D.12) 



(D.13) 



These two terms can be expressed in terms of magnetic and electric field gradients. 
Using equation (l90i) and (pp) from (I102p . we find 

(D.14) 



R 



b x R 



90 



2/ '(b x v R b x b) : (v R b) T - ^|v R± 5| 2 , 



B 



2B 2 < 



where M is the transpose of the matrix M. Employing that b x VRb x b 
(V R± b) T - (V R • b)(I -bb), we finally find 

,2 



-.(1) 
R 



9r R } 

b x 

09 



U 2 fl 

~B 



V R b : V R b 



"^(VR-bf-^lVa^l 2 . (D.15) 



B 



2B 2< 
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Substituting this expression and the definition of Sp ^ in (I108p into equation (1D.13|) gives 

,2 



^-V R b : V R b- ^(V R b) 2 



// 



2B 



A 2 



2B 
Vr 



Vr,S| 



2X 2 B 2 \ ( R±//Ae ) 



45 2 
b x V (R± /A. 



Using equations ( 147|) and ( 148|) . the term 72 in ( 1D.6j) becomes 

W^4 2) 



To = b ■ r 



R 



dudO 



0. 



(D.16) 



(D.17) 



Employing equations fj4^|) and fl5T?|) the term 73 in (ID. 71) is written as 




(D.18) 



Finally, the term 71 in flD.8|) is manipulated to yield 



Ta = — ( MMi 



Using equation (loUl) to write — uh ■ + BY 



+ BT { g ] ) and equation (SHD to write b ■ T 
(u\) = and (fii) = 0, we obtain 



(i) 

R 



T4 = ^ 



u d(uf) 



2 du V dfidO 

Using equations (T4"9|) and fl50|) . we find 



90 2 



-B{dSf/d6) + - «b ■ 
b • r^, and employing that 



(D.20) 



^1) p/ r (D g2 4 2) \ 




Here we have used 



a4 2) a 2 4 2) 



94 2) 9 2 4 2) 



09 dfidO du d9 2 



} dsf d 2 sf 

' { 89 8/286 




(D.21) 



(D.22) 



where we have integrated by parts in 9 to obtain the first equality. 

Substituting the results 7T« = A(<j>) + A0, f HXTB]) . (iDTfjl . flTXTSj) and rtD^Tj) into 
equation (ID.3[) gives 

2 «• ~ 3,u 2 



(2) 



2 

|^b • v R b • v R s + ^(i -bb) : v r v r b ■ b 



4fi 2 



VriB 



+ ^V R b : V R b 



ii 2 /i 
25 



(Vr ■ b) 



A 
A 



(R 2 - V 



(Rj./Ae)< 



25 

^^V R 5 • V(R ± /Ae)(0) - (V(R ± /Ae)$ • (b X V(R ± / Ae 

2\ 



(M 2 ) M (9(m 2 ) 

2 2 9m 



£? r 



(1) 9 2 4 2) \ Bd 
d^dO J 2d/i 




(D.23) 
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Employing equations (149[). (joTjl . (J52jl and f l69|) to write 

7(2) 



0^ [ 



00 



-A<f> - uui - 



(1) 



(D.24) 



we find 



p/ r (D^ 2 4 2) \ Bd 
\ 9 dfidO 1 2 8fi 



OS 



(2)' 



A 2 



<9(0 2 ) u 2 0{u\) 



Au I d(j) 
A 2 !? \ 0(n/\ 2 \ 



00 

Au 
~B 



2X 2 B0(/i/X 2 ) 2B Oil 



Oui 



-A 



or 



(i) 



or 



(i) 



dp, 



u ( U\- 



Substituting this result into (ID.23j) gives equation (I132p with and \f^£ as in (11331) 
and (11341) . and given by 

,2 



(D.25) 



,(2) 



(2) 
_0 



•|^b ■ v R b • v R 5 + -bb) : v r v r b • b 



4B 2 



\V R± B\ 



u 2 \i 



U 2 [X 



+ lT#VRb : V R b - -£(V R • b) 2 + 



25 



25 



(u 2 ) uQ(u\) u 2 0(u\) 
~2~ + 2~0u~ ~ 2B dfj, 



or 



(i) 



— U { U\- 



djj, 



(D.26) 



To obtain equation f II 3 5 j) from equation (1D.26j) we only need to calculate (u\ ) and 
(u\(Or^ /Ojj)). The gyroaverage of u\ is 



(«?> 



ii 2 /z 
IT 



R 2 

b.v R b| 2 + ^ 



(p x b) ■ V R b • p 



j(b-V R xbf, 



(D.27) 



where we have used the definition of u\ in f illip , we have taken the gyroaverage (pp) 
from (j!02p . and we have employed equation (jl07p to write [p(p x b) + (p x b)p] : VrI) = 
2(p x b) ■ Vr£> ■ p — 2(/i/£?)b ■ Vr x b and hence obtain 



|[p(pxb) + (pxb)p 

4/i 2 



: V R b 



4 



[pxb) - V R b • p 



B 2 



(b • Vr X b) 2 . 



(D.28) 



The gyroaverage of the second term in (1D.27P is given by 



[p x b) • V R b • p 



2B 2 



(b • Vr X b) 1 



25 2 



v R± b : (v R± b) n 



2 2 
/r (b x v R b x b) : v R b = -^-(b • v R x by 



2B 2 

+ ^V R ,b : (V R ,b) T - ^ 



2B 2 

(Vr • b) 2 , 



(D.29) 



where we have used 
,2 



(pipjpkpi) 



[(5ij - bibj)(8 k i - b k bi) + (8 ik - bib k )(5ji - bjk) 



P 

2B 2 

+(8u - bibi)(5 jk - bjb k )}. 



(D.30) 
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Here Sij is the Kronecker delta. We have employed b x VrIo x b = ( VR±b) T — ( Vr ■ b) ( I 
— bb) to obtain the second equality in (ID.29p . Substituting equation f)D.29|) into 
equation ( 1D.27j) gives 



(ul) 



u 2 /i 



b • V R b| 2 + ^V R± b : (V R± b) T - ^(V R • b) 

4 o 

-^(b-v R xb) 2 . 

^(i) 



The gyroaverage of ui(dT /<9/x) is 



Ux- 



^b • v R b ■ v R 5. 



(D.31) 



(D.32) 



Finally, substituting equations flD.31j) and (ID.32|) into f|D.26|) gives (I135p 



Appendix E. Poisson bracket 



In this Appendix we prove that the Poisson bracket that corresponds to the gyrokinetic 
Lagrangian in ( I136p is f )139p . Since the symplectic part0 of the gyrokinetic Lagrangian 



(I136p is exactly the same as in [25], the Poisson bracket in gyrokinetic coordinates will 
also be. 



As explained in subsection 13. 1\ to obtain the Poisson bracket, given in (|22|) . we 
need to calculate the inverse of the matrix L in (|20|) . We explicitly write this matrix 
by writing the gyrokinetic coordinates as {Z a } & a=1) with (Z 1 , Z 2 , Z 3 ) = R, Z 4 = u, 
Z 5 = fi and Z 6 = 9. The gyrokinetic Lagrangian f)136p is written as in f )29|) . with 
(r\,r 2 , T 3 ) = e^ 1 A(R) + ub(R) - e/iK(R), f 4 = 0, f 5 = and f 6 = -e/x. Then, using 
( 120|) for the gyrokinetic Lagrangian, we find that the matrix L is given by 

dT dT a 



L a /3 

or in matrix form 



oz a dz? ' 



(E.l) 



-e- i B*x I 



b , eK 



V 



b 

-eK 

6" 



' ' 

J- I 

' ' 

4- 1 

i i e 







— e 

"o 7 



(E.2) 



jj It is common to call the symplectic part of a phase-space Lagrangian to the piece linear in the time 
derivatives of the coordinates. 
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where B* is defined in (11401) . Its inverse is given by 

/•. ; : i 



B* 



ebx I 



cf 



B* 

~o 



V 



-eKxb . . . i — B* • K i —e~B 





0~~" 
0~~" 

-1 R* 



\ 



eKxb 



B* K 

e 



1 5* 



II. 







(E.3) 



/ 



where B* is defined in ( Il4ip . It is easy to check by direct calculation that P is the 
inverse of L. 

The Poisson bracket of two functions F(Z) and G(Z) is then given by equation (1221) 
that can be compactly rewritten as ( I139p . 



Appendix F. Calculation of the Jacobian 

In this Appendix we show that the determinant of the Jacobian matrix of the gyrokinetic 
transformation is B*, defined in ( 1141 p . This result coincides with the results in [25] 
because of our choice for the final form of the Lagrangian 



To obtain the Jacobian of the gyrokinetic transformation, we use the matrix L, 
defined in f[2"Up . This matrix is defined for both the original coordinates {r, v} and the 
new gyrokinetic coordinates {R, u,fi,9}. The matrix L in the original phase space and 
the matrix L in the new gyrokinetic phase space are related by the Jacobian matrix 
of the gyrokinetic transformation T e . It is possible to use this relation to calculate the 
determinant of the Jacobian matrix by calculating the matrix L in both coordinate 
systems. 

We denote the original coordinates by {X a } & a=l) with (A 1 , A 2 , A 3 ) = r and 
(A 4 , A 5 , A 6 ) = v, and the gyrokinetic coordinates by {Z a } b a=1 , with (Z 1 , Z 2 , Z 3 ) = R, 
Z 4 = u, Z 5 = fi and Z % = 9. The Jacobian matrix of the transformation is given by 

dX a (Z) 



)?(Z) 



dZ? 



We write the Lagrangian in the coordinates X as 

6 J VOL 

£ x = £^(x)^-# x (x,t), 



a=l 



dt 



(F.1) 



(F.2) 



with (71, 72, 73) = e _1 A(r) + v, 74 = 0, 75 = and 76 = 0. We write the Lagrangian C 
in gyrokinetic coordinates Z as in f[2"9"p . with (r 1; T 2 , T 3 ) = e _1 A(R) +ub(R) — e/iK(R), 
T 4 = 0, T 5 = and T 6 = — e/i. 

From the Lagrangians £ x and C, we obtain the matrix L, defined in (!2Up . in both 
coordinate systems, given by flE.ip for the gyrokinetic coordinates, and by 
, 07/3 <97 Q 



af3 



dX a dXi 3 



(F.3) 
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for the original coordinates. It is immediate to check that L = ££lJT e , with the 
superscript T standing for matrix transposition. It is then obvious that 



det(J T .) 



; det(L) 
det(Z) : 



(F.4) 



where we have used that the Jacobian of T t is positive at e = to determine the sign 
in front of the square root. Then, to calculate the Jacobian is enough to calculate the 
determinants of the matrices L and I. 
The matrix I is 



(F.5) 



-e _1 Bx 






I I 


- I 








I 








and the matrix L was given in flE.2j) . The determinant of I is 

det(Z) = 1. 
As for L, given in (IE.2I) . we have that 



(F.6) 



det(L) = e 2 



-e _1 B*x I 



b 







(F.7) 



Writing the matrix in the reference system {ei, £2, b}, where b = (0, 0, 1), we find that 

e" 1 ^ 



det(L) = e 2 






e _1 B* b 


-e^B* • e 2 





-e _1 B*-b 





e _1 B*-ei ! 







-e _1 B* ■ ei 


! 


-1 








i 

1 1 






~ l B* 



{B 



*\2 



Substituting this result and flF.6}) into (IF.4|) . we finally obtain 
det(J T J = B\ y 

Appendix G. Conservation of phase-space volume 



(F.8) 
(F.9) 



In this Appendix we prove that 



E 

a=l 



d 



det(J Te ) 



dZ° 
~dt 



(G.l) 
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This equation is satisfied by any gyrokinetic Lagrangian C with a symplectic part r Q 
that is independent of time, as is in our Lagrangian f)136p . Relation (10. lj) gives equation 
([153]). 

To prove f lGTTj) we use equations (JSTj), (Q) and (JEI]). From (ED) and (El we 

find 



E 

a=l 



g 



det(J Te ) 



dZ° 
~~dt 



1 



2^detW v 
+2det(L)]T ^ 



dt «9Z Q 



dt 



(G.2) 



and using f[2Tj) for the gyrokinetic Lagrangian £, we obtain 

<9# 



dZ° 
~dt 



a = 1, 2, . . . , 6. 



0=1 

Employing that (L~ 1 ) al3 = — (L~ 1 )i 3a , we find that this equation leads to 
d fdZ a \ _ A diL- 1 )* 13 dH 



(G.3) 



E 



<9Z Q V dt 



(G.4) 



a=l x / a,/3=l 

Equation flC4|) can be further simplified by using (IA.4j) . giving 

6 



E 



d fdZ 



dt 



^-i^ap dLpj dZ 1 



a,/3,7=l 
6 



<9Z Q (it 



s E(^)""(§? + 

a,/3, 7 =l V 



dL^ a \ dZ 1 



dZ? / dt 



(G.5) 



where to obtain the last equality we have used that L a p = —Lp a and (L _1 ) Q/3 = 
— (L~ 1 ) ,3q . Substituting ( 1G.5I) into flG.2|) and using that the derivatives of a determinant 
are 



VetOL) = det(L) £ (L- 1 ' 

(8,7=1 



(G.6) 



finally gives (IG.lj) . To prove that all the terms cancel we just need to use that, as 
trivially deduced from its definition ( f2"U|) . L a/3 satisfies ( lA.ip with n = 3. 



Appendix H. Proof of equation (1184ft 

In this Appendix we prove ( I184p by showing that det( Jz p0 ,_>z p (Zpo)) = 
By p (Z p0 )/B^ p (Zp(Z p o, to; t)). First we evaluate the time derivative of det( Jz p0 m.z p ), 
given by 

d 6 6 dZ 13 

-[det(J Zp0 ^z p )] = det(J Zp0 ^ Zp ) E ( J z^z p )l^Jr = det(J Zp0 ^z p ) E 

/3, 7 =1 P° 0=1 
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where we have employed the formula for the derivative of a determinant in flG.61) and 
omitted the arguments hoping that no confusion will be caused. Using (I153p . or its 
equivalent (IG.lj) . we obtain that 



^—^(J^)] = ^.^.^ IWW))!, (H.2) 

that is, the product B*^ p (Z p (Z p0 , t ; t))det( Jz p0 ^z p ) is constant in time. Since for t = t 
the map Z p (Z p0 , t ; t) is the identity, giving det( Jz p0 n-z p ) = 1, we find that the constant 
is Bl p (Z p0 ), implying that det( J ZpO ^z p (Z p0 )) = ^(Z^)/^ Z p (Z p0 , t ; *))• 

Appendix I. Manipulations leading to equation (1186ft 

In this Appendix we obtain equation (I186p from ( I185p . 

First, we evaluate the variations with respect to y?(r,t) in dl85[) term by term. For 
the first term in (11851) . we find 

J dtd 3 r V5ip(r,t) ■ V(p{r,t) = -J dt d 3 r 5ip(r, t) VVO, t), (1.1) 

where we have integrated by parts and we have taken into account that Sip = at the 
boundary. Using relation (ITUj) and equation (I162p to write 

4> P (R P , n P , 9 P , t) = J d 3 r 5 (r p + j-p(R p , n P , Op) - r^j <p(r, t), (1.2) 



we find 



5 v (0 p (Rp, fMp, 9 P , t)) = J d 3 r 5ip(r, t) (^5 (r p + j-p ~ r ^ 



(1.3) 



Using relation ( II. 2 p again and the identity 

^^(V (Rp±As) $ p • (b x V (Rp±As) p )) = ^(V (Rp±As) (c^$ p ) ■ (b x V {Rp±As) p )> 
-^(V (Rp±As) (^0 p ) • (b x V( Rp±As) i p )) 

= -(V (Rp±As) (^0 p ) • (b. x V (Rp±As) $ p )>, (1.4) 

where to obtain the last equality we have integrated by parts in 9 P and we have used 
that d^ p /dOp = (ftp and that d(5 ip Qp)/d9 p = b^p, we obtain 

p, Up, [lp, 9 P , t) 

= J d 3 r5v(r,t)(^-j-V Ilp (Rp,Up,fip,9p,t) ■ V (RpxAs) 5 (r p + j-p - r^j 

d 



+ ^/ip (Rp; Up, lip, 9 p , t) - 



5(R p + f p-r 

A p 



d/jp 

+V 9p (Rp, u p , Hp, Bp, t)5 (r p + yp - r^j 



(1.5) 



In this equation we have separated the different terms into three types: the terms that 
contain the gradient of the delta function, the terms that contain the derivative with 
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respect to \x v of the delta function, and the terms that contain the delta function. The 
coefficient multiplying the gradient of the delta function is 



V Rp (R p , Up, n P , 9 P , t) 



A u " 

" b x V iKpM % - x V Rp b ■ p 



\B 2 



2B 2 Il p ± 



- — [pp-(pxb)(pxb))-Vn p B; 



(1.6) 



the coefficient multiplying the derivative with respect to \i p of the delta function is 



V^ p (Rp,Up,//p,e p ,t) = --g^p - -gb • V Rp b • p 

U r , 



+jV Rp b: (p(pxb) + (pxb)pj; 
and the coefficient multiplying the delta function is 



(1.7) 



Ve p (&p, u y /V 



A, 



VrB • P 



V Rp b: p(p x b) + (p x b)p . 



XlB d^p/Xl) B 



2fx p B 



b • V Rp b ■ p 



(L* 



Substituting CI) . (Ojl and (JH5J) into f lT85|) . and employing (1X83]) and f lTMj) to write 
the integrals as integrals over Z p and not over Z p0 , we find that the variations are of the 
form given in (I186p with 



V(r,t) 



X 2 c 

A De t s -r-,2 



L 2 



W(r,t) 



^Zp J d 3 R p dupdn P d9p B^ p (R p , u p , /i p )F p (R p , u p , /i p , t) j<5 ^R p + j-P ~ r 



+ 



V R „ ■ V R J r, 



A p 



5 I R p + — p - r 

Ar, 



+ 5 R 



A, 



(1.9) 



Note that we did not need to keep the gyroaverages in (II. 3j) and (JL5J) because there is 
an overall integral in 6 P and neither B* nor F p depend on 9 P . 

We have only left to prove that ( II. 9 P is equal to (j!87j) . Using the result in ( II lOj) 
and employing ( I102p and ( I107p to write 



-g-V Rp±J B + - 



pp- (p x b)(p x b) • V Rp 5 = pp • V Rp S, 



(1.10) 



it is clear that V Rp is exactly the perpendicular component of R p> 2- Using (I112p . we 
find that U Mp = It only remains to use that Vq p = —d9 Pt i/ 'd6 p , with 9 Pj i given in 
( I113p . To obtain this identity we have used p = d{p x h)/d9 p , p{p x b) + (p x b)p = 
-(l/2)d[pp - (p x b)(p x h)]/d9 p and P = d%/89 p . Using U 6p = -89 p>1 /d9 p and 
integrating by parts in 9 p in ( II. 9 P gives the final form ( I187p . 
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